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NUMERICAL SOLUTION OP THE COMPRESSIBLE- LAMINAE, TRANSITIONAL, 

AND TURBULENT BOUNDARY LAYER EQUATIONS WITH 

■COMPARISONS TO EXPERIMENTAL. DATA - - 

By 

Julius Elmore Harris 
ABSTRACT 

A system of eq.uatioii.s describing the laminar, transitional,, and 

* • * ' ^ 

■ turbiolent compressible boundaiy layers for either planar- or axi symmetric 
flows together with a numerical method by .which the > system can be 
accurately solved is presented. Stability theoiy and its relation to 
transition is discussed, and methods are presented with which reasonable 
estimates may be made of the location of transition and the extent of 
the transitional flow region. These methods are used in the numerical . 
procedure to connect the three separate boundary layer regimes into one 
continuous region described by one system of .governing equations. The 
transitional boundary layer structure is developed from a consideration 
of the statistical production of turbulent spots. ' The fully developed 
turbulent region is treated by, replacing the Reynolds, stress terms with 
an eddy viscosity model. A specifiable turbulent Prandtl number is used 
to relate the turbulent flux of heat to the eddy viscosity. 

The momerical method used to solve the system of equations is a 
three-point implicit finite -difference scheme for variable grid-point 
spacing in both spatial coordinates. The method is self starting; that 
ia, it requires no experimental data input, and is highly efficient, with 



regards to flexibility.^ computer processing time, and accuracy. The • 
method is inherently stable; no constraint is placed on the system by 
a grid-point spacing stability parameter. To the author's knowledge 
this represents the first time this paorbicular numerical procedure has 
been applied to transitional and fully tiirbulent boundary layer flows 
■as well as the first time the transitional flow structure has been 
included in such a procedure. 

Numerous test cases are presented and the results are - compared' with 
experimental data for supersonic and hypersonic flows. These test cases 
include flows with both favorable and adverse pressure gradient historie 
mass flux at the wall, and transverse curvature. The results clearly 
indicate that the system of equations and the numerical procedure by 
which they are solved can be used to accurately predict the character- 
istics of laminar, transitional, and turbulent conpre s sib le boundary- 


layer flows. 
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VI. UJTRODUCTION 


The 'boundary layer concept first introduced by Prandtl (ref. I) 
in 190^ divides the flow field over an arbitrary surface into two 
distinct regions; an inviscid outer region in which solutions to the 
Euler equations describe the flow field characteristic and a viscous 
inner region where the classical boundary-layer equations are valid. 

The boundary-layer region may be further divided into three categories; 

■ ' i 

namely, laminar, transitional, and turbulent. 

The laminar boimdary layer has received considerable attention 

ji ' f - 

over the past 60 years, especially over the past decade with the 

increased availalaility of high-speed digital computers. Early solutions 

of the boundary layer equations were obtained for incompressible, zero 

pressure gradient, adiabatic flows These solutions' were generally 

obtained by series expansion techniques and -were elegant in formulation, 

* > 

but time-consuming to obtain; however, they did yield valuable information 
for extrapolation to more complex flow systems. The addition of pressure 
gradients into the problem required the use of nunerical techniques which 
were impractical without high-speed digital computer systems. This 
difficulty led to the development of momentum integral approaches in 
which the actual boundary layer equations were not satisfied at each 
point across the viscous region, but instead were .satisfied only in the 
average. Reviews of these early methods are given in references' 2 , 5 ; 
and i*-. 


1 
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As flight speeds increased, it became' necessary to include the 
effects of compressibility. The inclusion of compressibility coupled 
the momentum and energy equations and resulted in a formidable system 
of five equations, three of which .(the conservation of mass, momentum, 
and energy) were nonlinear partial differential equations. The 
requirement of simultaneous solution of this system prevented any 
complete numerical solutions from being obtained until the advent of 

j ' - 

early digital computer- systems. There then appeared a number of similar 
and so-called local similarity solutions, i A - review of these methods 
and solutions is presented in reference 5 * Finally^ in the early part 
of the past decade the complete , nonsimilar laminar equations for the 
compressible, nonadiabatic boundary layer were solved to a high degree 
of accuracy by finite difference techniques (see Blottner, ref. 6). 

One of the earliest observations made by students of fluid mechanics 
was that, in general, a turbulent or random-like state of motion was- the 
most natura-l state of fluid flow. 0. Reynolds (ref. 7 ) in his now 
classical study of pipe flow observed that at some value of the parameter 
Ud/v, where TJ, d, and v represent the mean velocity of the flow, 
the diameter of the pipe,, and the kinematic viscosity, respectively, 
the laminar flow degenerated to a turbulent state in which the molecular 
viscosity became of secondary importance in relation to the Iclnematics 
of the flow. Early investigations of turbulent flow categorically 
divided nonlaminar flow into two regions; transitional and turbulent. 

Transition, and the ability to accurately predict its location 
on an arbitrary surface has been the object of intensive analytical 
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and experimental research, for many years. .A complete imder standing of 
the transition process as well as the ability to predict its location 
for general flow configurations has not yet' been achieved. However^ 
the effects of basic flow parameters on transition have been studied 
in detail. The flow within the transition region, which is basically 
nonlinear in character is neither fully laminar nor fully turbulent but 
a combination of both. The detailed mean flow within the transition 
region itself has not. been studied as extensively as the location of 
transition and the characteristics of the eventual fully developed 
turbulent' boundary layer. Consequently, little if any effort has been 
specifically directed towards the problem of developing a suitable 
system' of equations that would describe the mean characteristics of 
transitional flow. There have been a few ejcperiments in which the mean 
profiles were . measured as well as some where the inteimiittent character 
of the flow was studied, but more detailed work is still required. 

There 'have, of course, been many experimental tests at high speeds ,in 
which the heat transfer at the wall has been measured, but this describes 
little if anything of the flow structure away from the wall. Savulescu 
(ref. 8) has recently presented one of the first thorough reviews of 
transition phenomena. Mofkovin (ref. 9) recently completed the most 
current and extensive 'reyiew of modem stability theory and experimental 
hypersonic transition. The characteristics of transitional boundary • 

i 1 • I ' 

layers for low speed incompressible flow as' well as compressible flow 

' ' • ' f '> 

(ref. 10) has; received some attention. 'These results at least allow 
■workable models of the me^ flow structure in the transition region to 
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be formulated and applied tentatively to compressible flow systems. 

It appears that at the present time it is not possible to obtain 
completely general solutions for transitional flow. However ^ it is ' 
possible to develop models for the mean flow from existing data on the 
intermittency distributions which represents in a broad sense the 
statistical distribution of tirrbulence intensities. 

Compressible turbulent boundary-layer flows have received much 
attention over the past decade because of the vital need of being able 
to accurately predict heat transfer and skiu friction for hi^ performance' • 
aerospace vehicles. However, most of the work has been experimental with 
the main objective directed towards developing empirical or semiempirical 
correlation techniques. Little if any effort was devoted to obtaining 
numerical solutions of the 'equations for turbulent boundary layers until 
a' few years ago. The principal difficulties were associated with the 
modeling of the turbulent transport terms as well as techniques for 
obtaining solutions on existing, digital computer systems. Even today 

f, * , ♦ 

because of the limited understanding of these turbulent transport 

processes, completely- general solutions of the mean turbulent boundary 
; \ 

layer equations are not possible . However, by modeling the turbulent 
transport terms through eddy viscosity or mixing length concepts it is 
possible to solve the system of equations directly. Reviews of recent 
analytical advances are contained in references 11 and 12 for 
incompressihle flow and references 13 and l4 for compressible flows. 

The purpose for the present thesis is to present a system of 
equations describing the laminar, transitional, and turbulent compressible 
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boundary layers and a solution technique with, wbich. the system may be 

accurately solved for either planar *or -axisymmetric perfect gas flows. 

The solution technique has "been 'found to yield accurate* results for 

compressible laminar^ i;ransitionai_, and fully developed turbulent 

boundary layers with pressure gradients^ heat transfer, and mass 
’ < * ’ 
transfer at the wall. The solution technique utilizes 5-Point implicit 

difference relations and the method first developed by Flugge-Lotz and 

Blottner (ref. 15) and later improved upon'by Davis and Flugge-Lotz 

(ref. l6). to solve the difference equations. The equations are solved 

in the transformed plane,. Transverse curvature terms are retained, 

and variable entropy effects may be included. The transition region is 

modeled by utilizing an intermittency distribution which describes the 

statistical distribution of turbulent spots and modifies the models of 

the turbulent transport processes.' The turbulent region is treated 

by solving the mean turbulent boundary layer equations in which the 

Reynolds stress terms are replaced by an eddy viscosity model, and a 

specifiable turbulent Prandtl number function- relates the turbulent 

flux of heat to the 'eddy viscosity. The eddy viscosity model is based 

upon existing experimental data. 



VII. MATHEMATICAL DESCRIPTION OF THE EftUATIOITS FOR THE LAMIHAR, 


TRANSITIONAL, AND ' TURBULENT COMPRESSIBLE' 

BOUNDARY LAYER 

^ f I 

In this- chapter the governing equations 'for the compressihle 

boundary layer togethe-r with the required boundary conditions are 

' - ** 

presented. Special attention .is devoted to the eddy viscosity and eddy 

‘ I ' 

conductivity models used to represent the apparent turbulent shear and 

heat flux terms .appearing in- the mean turbulent boundary-layer equations. 

Boundary layer .stability, transition,' and transitional region flow 

» 

structure are also 'discussed. 

t 

7.1 The System of" Par'bial Differential Equations 
T-l.l Geometry and Notation 

The orthogonal coordinate system chosen for the present analysis 
is presented in figure I. The boundary layer coordinate- system is 
denoted by x* and y* which are tangent to *and normal to the surface, 
respectively. 'The origin of both the boundary layer coordinate system, 
(.X' 5 <-.,.y*)' and the body coordinate .system,- is located at the 

stagnation point for blunt body flows as shown in figure, 1 , or at the 
leading edge for sharp-tipped cones or planar surfaces. The velocity 
components u* and v* are oriented in the x-^^- and y* direction, 
respectively. Transverse curvature terns are retained because of their 
importance in the development of boundary-layer flows over slender bodies 
of revolution where the boundary-layer thickness may become of the order 
of the body radius, r*. A discussion of transverse curvature and its 
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effect on the resulting' hoimdary layer characteristics is presented 
hy Lewis (ref- !?)• The radial coordinate, r* represents the true 
radial distance from the centerline of the hody to a given point (x*,y-*) 
in the boundary layer. The angle 0 is the angle between the z* axis 
and local tangent evaluated at (x*,o). The coordinates aiid 

(x* „^o)' represent the location at which transition is initiated and 
completed, respectively. 

7.1.2 Differential Equations 

The flow of a compressible, viscous, heat conducting fluid is 
mathematically described by the continuity, Navier- Stokes, and energy 
equations together with an equation of state, a -heat conductivity law, 
and a viscosity law. For flows at large Reynolds numbers, Prandtl 
(ref. l) has. shown that the Navier-rStokes and energy equations can be 
simplified to a form now recognized as the classical boundary layer 
equations. These equations may be written as follows (see ref. 2);. 

Continuity * ■ 

» *• * , 

*, - — (r.*'^p'%*) +■ (r*^p*v*) = 0 (7*1) 

dx^^ ' . ’ 3y* ■ . 

* ‘ F 

Momentum 



(7.2) 
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Figure. 1. 


Coordinate system and notation. 
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( 7 . 5 ) . 
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Energy 


1 * 


u* + (v* + 

Sx* P V p.* idy* P ! 


<ix* v*J ^y* 


&u* ' 


^ 0* Sy*^ P J \^ay* 


<)1 + 


2 -- 


+: 


1 5 . 

r*j 3y*‘ 


c|p*v*'T*'‘) 


- ■p*u*’v*’ 


Sy* 

( 7.60 


The mean turbulent equations are identical to those for the 

laminar boundary layer (eq,s. (T*!) to (7*3)) with the exception of the 

correlations of turbulent fluctuating quantities which, are the apparent 

■masS; shear j and heat flux terms caused by the action of turbulence. 

The main .problem of calculating tiirbulent flows from this set 

of -equations is concerned with how to relate these turbulent correlations 

to the mean flow and thereby obtain a closed system of equations . In 

the present analysis, the apparent mass flux term, the apparent 

shear stress term,, p-^ u*’v*' (Reynolds stress term) and the apparent 

heat flux term, C* p* v*'T*' are modeled or represented by a new 

velocity component, v*, an eddy viscosity e*, and an eddy conductivity 
* 

Kqi , re spe c t ively . 

A new velocity component normal to the surface is defined as 

follows: 

V* = V* + 

p'* 


( 7 . 7 ) 
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The eddy viscosity is defined as 


e* = - p* 


u*'v*' 


(7.8,) 


and the eddy conductivity as 


Kr^ = -- C* p* 




7T ■ -P r 


(7.9) 


The static tunMlent Praiidtl number is, defined as, follows:. 


>*'v*',/STV8y*) 


*^t 'i ■ 


v*'T*' \du*/5y*J 


(7.10) 


Equation (7.I0) can then be expressed in terms of, equations (7.8) and 

(7.9) as 


■* * 
■Cp€ 

^t =-pr 

% 


(7.11) 


The- eddy viscosity model used in the present analysis is discussed in 
Section 7.^'. 

In terms of equations (7*7) through (7.II) the governing 
differential equations may be written as follows: 

Continuity 


T^(r*^p^*) + :^(r*'^p-5«v*) == 0- 
ox* dy* 


(7.12) 
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Momentum 


p* + V* 

\ Sx* -By*/ 


^ H. 

<3X* 3 ,*j Sy*\ Sy*y* 


(T-15) 


Energy 


pX» ^C|T*) -h ^0|T*) 


. u* ^ 


Sx* 

_1 ^ 

r*j ^y* 


' \Sy*/ 






(7.14-) 


The terms e and e appearing in equations (7-13) and (7*1^) are 
defined as follows: 



(1 + ^ r) 


( 7 . 15 ) 


and ' 

I 

>> 


/V 

e 


/ e* cf ^ 

= (1 + ^ r 

p* • 


(7. 16 ) 


. { 


respectively. ‘ ^ 

* ■ ^ 

Tlie ’function, T, appearing in equations ( 7 * 15 ) and (7*l6) 

represents the, ,streanMise intermittency 'distribution in the transitional 
• ■ . . ^ 

region of the boundary layer, r assumes a value of zero in the laminar 
, « 

region of the boundary layer. and a value of unity in the fully turbulent 
boundary layer. The variation of r -within the transitional region of 




the houndary layer depends upon^the statistical growth and distribution 
of turbulent spots. The intermittency function is discussed in 
Section - ' ' 

In order to complete the system of equations the perfect gas law 
and Sutherlands viscosity relation are introduced. 

Gas Law 



■ (7.17) 


Viscosity Law 





(7.18) 


The system of governing euq.ations then consists of three nonlinear 
partial differential equations and two algebraic re'latxons. Two of 
these differential equations (eqs. (7.13)' and (7.14)) are second order 
while the remaining differential equation (7*12) is first order. 
Consequently, if suitable relations for e*, and P can be specified 

there are five tmlmowns.; .namely, u*, v*, p*, T*, and ji* and five 

equations . 

The pressure gradient term appearing in equations (7-i3)' aiid 
(7-14) is replaced by the Bernoulli relation; namely 


= - p*u* 

dx* ® ® dx* 


(7-19) 



•lit 

which is determined from an inviscid solution. If -variable entro^r 
is considered (see Section T*9)' dp*/dx* is retained in equations (7 •15) 
and ( 7 . lit-). 

7-1.3 Boundary Conditions 

In order to obtain a unique solution to the system of .governing 
equations it is necessary to satisfy the particular boundary conditions 
of the problem under consideration. These conditions are shown 
schematically in figure 2 . 

u"(x*,ye*) = Ue*(x*n 


■) 

)' 

^Figure 2.- Boundary conditions in’ the physical plane. 

‘ .V ’ ^ • - ' . • ' 

The velocity, and temperature distribution at the edge of the boundary 

• > -f V .■■■.* 

< ' t 1 ■ ; ■ ' 

layer afe determined from the shape of .the body by using inviscid flow 

t 

theory and are discussed in Section '7»8‘.'; The no- slip condition is- 
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imposed at the wall; however, arbitrary distributions of v* and 1* 
may be specified. 

The parabolic nature of equations (7-15) and (7»1^) requires 
that the initial velocity and temperature profiles be specified at x 
These initial profiles are obtained in the present investigation from 
either an exact numerical solution of the similar boundary layer 
equations or from experimental data and are discussed in Section 8.1. 5- 


7.2 The System of Non-Dimensional Equations 
7-2.1 Non-Dimensional Variables 

The non-dimensional variables chosen for the present analysis 
are defined as follows: 

Dependent Variables 


u = u*/u* 
ref 

p -.pVpj^f 

T = T*/t* 
ref 

u. = u*/ u* 

^ ^^ref ■ 


Independent Variables 


X = xf/L* 

y = y*/L^ ) 


(7. 20a) 


r 


rVL^J 


(7. 20b) 


H-ifc 


l6 


The reference values -of density and velocity used in equations (7- 20) 
are tafeen to be those of the free stream^ the reference temperature 
is taken to be and the reference viscosity is the value of 

the jnolecular viscosity evaluated from the Sutherland law at the 
reference temperature. The reference length (characteristic length) 

L* may be taken as the nose radius, for example, in the case of a 
spherically blimted body or as any unit length for planar flows. (See 
equations (,7*125.) bo (7-128).) 

7*2.2 Non-Dimensional Governing Equations 

’• ■ f , * 

The non-dimensional equations are obtained by substituting 
•equations (7*20) into equations (7.12)-, (7 -13), ' and (7*1^) and are as_ 
follows ; 


Continuity 


::^(r'^pu‘) + t;^r'^pv) = 0 ' , 

ox . oy 


(7.21) 


Momentum 






Rref 


1 87 d — dm 


( 7 - 22 ) 


Energy 


/ St , - 3t\ dp 1 

p u — + V — I = u ^ + — ~ 
\ dx, dy/ dx 




_1 _d 

dy\” 0 dyj 


+ ^ — .frJ HI 


(7*25) 


The parameter R^-gf appearing in equations (7*2l), (7*22), and (7-23) 
is a reference Reynolds number defined as 
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Pref^ref^ 

Rref = 

IJ^ef 


( 7 . 24 ) 


7.2.5 The Stretched Coordinate System 

In order ho remove the reference Reynolds number from hhe governing 

* f 

non-dimensional’ equations as well as stretch the coordinate normal to 
the wall a new independent variahle^ Y, is defined; that is 


Y = t/o>' 


(7.25) 


A new velocity, component is ^sO' defined by the relation 


. . ■ + v 

. V = — 

^ . cn 

* i 

where the parameter cn is defined as follows: 


(7.26) 


“ (^ef^ 


(7.27) 


The non-dimensional governing equations then become in terms of 
the stretched variables 


Continuity 


Momentum 


rr^r^pu) + ;^(r'5 pv*^) = 0 
ox oY 


(7.28) 


/ 5u + 5u\ dp 1 d / -i _ 8u' 
V dx dY,/ dx- rJ ,dY\ dY 


( 7 ^ 29 ) 


Energy 


/ dr + dr' 
p u — + v+ ~ 
\ dx dY 


dp 1 d / ^ p€ dr\ , — /du\ 

dx rJ dY\ P dY/ IdY/ 


(7. 30) 
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The equation of state (eq. (7*1^7)) and the viscosity law (eq. (7«l8)) 
may he expressed as follows; 

Equation of State 



Viscosity Law 


(7.51) 





1 + S 
T '+ S 


(air only) 


(7.32) 


where S = 

The houndary conditions for the system of equations in the 
stretched^ non-dimensional plane are presented in figure 5* 



Figure 3*- Boundary conditions in the stretched plane.’ 
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7. 3 The System 'Of Transformed Equations 
The system of governing equations is singular at x = 0. The 
Erobstein-Elliot (ref. 20) and Levy-Lees (ref. 2l). transformation can he 
utilized to remove this singularity as well as control the. growth of 
the boundary layer as the solution proceeds downstream. In reference l6,, 
Davis and Fliigge-Lotz obtained solutions to the laminar botmdary layer 
equations in the real plane for axisymmetric flows. However, the 
solutions were obtained for only a few nose radii downstream of the 
stagnation point and the boundary layer growth problem was not serious. 

If solutions are required many nose radii downstream, say on the order 
of 1,000, then the boundary layer growth can cause problems if the 
solution is obtained in the physical plane unless provisions are made 
for checking and adding additional grid-point when required. 

The transformation' utilized can be- written as follows: 
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The relation between derivatives in the old (x>Y) and new 
coordinate systems are as follows; 


3\ -2i/ ( Sri\/ S\ 

/y \ Vti \ /\ 


a\ . W^/IY-S'l 


'x V5T 


(T.35) 


Two new parameters,^ ' F and 0 are introduced and are defined as 


F = u/u, 


■0 == t/t, 


(T.56-) 


as well as a transformed normal velocity 


2 1 /Sti\ pv'^r itj I 

V = ^ 

pup r^j V^j X/si - 

*^6 e'^e o ' ' 


(7-5T) 


The governing equations in the transformed plane can he expressed 


as follows: 


■Continuity 


3v BF 

— +2|~ + F= 0 

Bn, 


(7.,38.) 


Momentum 


Bf Bf B/, 2i,-r Bp\ , 2 ^ 

2|F — + V (t "’le — + p(F 0^ = 0 

B| Bt] BqV Bti/ 


(7-59) 
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Energy 


0| OT) OT^ \ O’ dT\j 


\ori j 


( 7 .^) 


where 


2 = (:pu)/(pu.)g 

a ^ u|/Tg 






CT.41) 


The parameter 2 can he written by using' the viscosity relation 
(eq. (7-52)) and the equation of state (eq. ,(7*31)) as- 


2 



1 + S 

L© + S 


(air only) 


( 7 .^ 2 ) 


where S •= S /t© • 


The transverse curvature term can be , expressed in terms of the 
transformed variable s as ’ , * ’ ’ 


t = ± 




(7.i)-3) 


where the t sign in. equation ■(7*^3-)' is required in order to- obtain 

« t 

t , 

solutions for axi symmetric boundary layer flows p'ver a body or inside 
a duct (nozzle), respectively. 
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The physical coordinate normal to the wall in the stretched, real 
plane is ohtained. from the inverse transformation; namely 


/l + 2 ViT" ti) cos 0 r'^ '\ 

- 1 i / Q 

V \J >v J 


Pe^e^o 


1/2"]; 


(7. 44) 


The selection of the correct sign in. equation (7*44) is made on the hasis 
of + for axisymmetrlc flow over bodies of revolution and - for flow 
inside of axisymmetric ducts (nozzles) . / ' ' 

The y- coordinate in the physical plane is obtained from 
equations ( 7 . 25 )^ ( 7 . 27 ).^, and (7-44); that is, y = coY. 

The boundary conditions in the transformed plane are as follows: 
Wall Boundary 


F(|,0) - 0 

V(g,0) =V^(|) 

0(1.0) =0j|) 


(7. 45a) 


■Edge Condit ions 


FU.Tie) =1 

©(e.Ti^y = 1 


(7.45b) 


The boundary condition at the wall for the transformed V component can 
be related to the physical plane as 


^w = 


_ J^/Pw'^w 


cod^r J \Pe^e 
e o 


(7.46) 
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where the no-slip constraint has heen imposed on .equation (7*5T).* It 

' » i 

should he noted that the apparent mass flux term appearing in 
equations (7-^) and (7-5) is zero at the wall. Therefore equation (7-^) 
can "be expressed in terms of the physical mass flux at the wall as 


V, 


w 


_ \/^ / Pw^w \ 


(7-^7) 


7 .'^ Eddy Viscosity and Eddy Conductivity 'Models 
7 A.I Physical Plane 

The turbulent boundary layer can be treated as a composite layer 
•consisting of an inner and outer region as shown schematically in 
figure 4.- (see Bradshaw, reference 22.) 
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The thickness of the Innei' region ranges between 0.15 to 0.26. 

The character of the flow depends primarily on the shear stress at the 
wall^ and the molecular viscosity, p.. The mean velocity distri- 

bution in this region responds quickly to changes in the wall conditions 
since the turbulent eddy structure is small. The mean velocity 
distribution is characteristically described by the "law of the wall" 
relation 





The law of the wall relation’ was first derived, by Prandti (ref; 25') . 

If the wall is smooth the inner region will .’contain a sublayer., usually 

referred to in the literature as, either the ’laminar sublayer or the 

viscous sublayer, adjacent to the wall. / , ■ ' _ - 

> 

The viscous sublayer is very thin in relatiqn to the total 

boundaiy-layer thickness. -The thickness of the sublayer usually ranges 

be-tween 0.0016 to 0.016. The layer is highly viscous in character) 

■consequently, the mean velocity is a linear function of* y. For the 

viscous sublayer, equation (7.48) can be written as 

* * 


u 


y Hr 




V 


* 


(7.49) 


The relative thinness of the viscous sublayer, and its importance in 
determining the correct shear stress and heat flux at the wall requires 
that very small grid-point spacing in the y-direction be utilized in 
the' wall region for the numerical procedure used in the present 
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analysis. This particularly important point is discussed in 
Section 8.1. 

The outer region of the turbulent bo-undary layer is very thick 
in relation to the wall region. It ranges in thickness from 0.86 to 0.95* 
The flow in ttiis region is basically independent of the molecular 
viscosity, dependent upon the wall shear stress*, and strongly affected 
by boundary layer edge conditions such as pressure gradients in the 
streamwise direction, dP/dx. The mean velocity distribution for the 
outer region is usually discribed by the velocity-defect law-: 


The basic argument leading to the, form of the velocity-defect law is 

t ■ ' ' ' 

that the reduction in velocity, -Ug<- ,u at a distance y . from the wall 
is the result of a tangential stress .at the wail, independent of how the 
stress arises but dependent on the distance to which the effect has 
diffused from the wall. The flow in the outer region shows similarity 
to wake flow, and the eddy structure is large scale in relation to the 
inner region. Consequently, the response of the mean velocity distri- 
bution to changes in boundary conditions is quite slow. The eddy 
viscosity model, then, must characterizes these basic flow structures if 
it is to yield accurate results. For additional discussion on the inner 
and outer regions see references 2 and 2h. 





Inner Region Model 


Tiie eddy viscosity model -used in the present analysis for the 

! 1 , • . 

inner region is based on the mixing-length hypothesis as developed by 
Prandtl (ref. 25) in I925. The eddy viscosity for this region 
referenced to the molecular viscosity may be expressed as follows: 




where 2.* the mixing lengthy may be written as 


(7.51) 


I* = K^y* ' (7.52) 


The value of Kp has been obtained experimentally and has a value of 
approximately 0.4^ the value which will be used in the present analysis. 
However^ Van Driest (ref. 26) concluded from an analysis, based upon 
experhnental data and the second problem of Stokes (ref. 27) (sometimes 
referred to as the 'Rayleigh problem^ ref. 28) that the correct, form for 
the mixing length in the viscous sublayer should be as follows: 


I* = 




1 - exp(- Jj) 


(T.53) 


where the exponential term is due to the damping effect of the wall on 
the turbulent fluctuations. The parameter A* is usually referred to 
as the damping constant.. The exponential term approaches zero at the 
outer edge of the viscous sublayer so that the law of the wall region 
equation, as expressed in equation (7-52) is valid. The damping 
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constant A* is a strong function of the wail houndary , conditions and 

Equation (7*.5^) was originally obtained for incompressible ^ zero 
pressure gradient^ solid wall flow,j that is _ p* = constant^ dp*/dx* = 0 
and = 0. The relation has^ however^ been applied to compressible 
flows where p* and v* are evaluated locally across the sublayer^ , 

A 

where p* is evaluated locally and an average value of v* is utilized, 
and where v* and p* are both evaluated at the wall conditions. In 
the pie sent analysis the mixing length will he defined as follows 


is defined as. 


A* = 26v4^\ 


\P^ 



( 7 - 55 ) 


_* 

where ' v is the average value of the kinematic viscosity taken over 
the viscous sublayer. The density and viscosity appearing in 
equation (7-5^) will be evaluated locally. The eddy viscosity for the 
inner region referenced to the molecular viscosity can then be written 
as 




[77^ *v 

2 

* 

rw y \ 



1 V* A*} 


Sy* 


{■7.56) 


where A* 


is defined in equation (7>5^)- 
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Cel^eci (ref. 29) has recently attempted to account for the effects 
of both pressure gradient and mass injection or removal at the wall on 
the damping constant . Cebeci begins with the two-dimensional x-momentum 
equation written as 


^ 8u* ^ Su* 
u* T — + V* — r— 
Sx* oy* 


p* dx* p* Sy*\ 


p*u*'v*'| 

8y* / 


(7.57) 


Cebeci then neglects the Reynolds stress term and assumes that in the 
wall region equation (?.57) can be written as (v* = O) 


■' . ■ '/ (7.. 58) 

dy* dx* * . 



It should be noted that for v* == 0 it is not necessary to neglect the 
Reynolds stress termj that is, r* in' equation (7*59) could be replaced 
with Tm. -However, if v* = 0 the Reynolds stress term must be 
neglected in the approach used by Cebeci. The wall shear stress term 
r* appearing in equation (7.5^) is then replaced by the local shear 
such that a new definition for the damping constant is obtained^ that is 

A* = 26 v.(-s- + sn ' 

yp* dx* p*y 


(7.60) 
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This approach is 'SimiTar to that 'followed by Patankar and Spalding 

i ^ \ . 

(ref. 50 ) where they utilized the Io.c^ .value' of, t* in equation (7*5^) 
instead of the' wall value ^ ‘ 

The assumptions made by Ceheci in obtaining .equation (,7*6o) are 
not necessarily valid. Furthermore j as previously 'mentioned^ the flow 
in the outer region of the turbulent boundary layer is predominately 
affected by pressure gradients and not the flow in the law of the wall 
region. Consequently, the pressure gradient correction to the damping 
constant as obtained by Ceheci is of questionable value and will not be 
utilized in the .present analysis. 

The effect of mass flux at the wall has also been treated by 
Bushnell and Beckwith (ref. 31 ) and Cebeci (ref. 29). Cebeci, for non- 
zero V* expresses equation (7*57) s-s 


dr* 

dy* 


~3(- 

V 

w 

T* - 

V* 


ap* 

dx* 


= 0 


(7.61) 


which can be integrated to- yield 


T* = 



+ 


dp* 


Vw 


dx*i 



(7.62) 


Equation (7*5^) can then be written as (where t* has been replaced 
with T*) 




A* = 26vV^ exp(^ y*| 
IP* \v* / 


i^dp*^ 

V* p*\dx*/ 




exp^^. y J 



( 7 - 63 ) 
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Buslmell and Beckwith (ref. 51) studied the effect of wall mass 
flux on the damping constant from experimental data for zero pressure- 
gradient flow. These results agree well with equation (7-63) for zero 
pressure-gradient flows; that is 

-1/2 

( 7 . 61 ^) 

Two apparently conflicting ideas by the author of the present 
paper should now be discussed concerning the development of Cebeci 
(ref., .29). Firsts it has been stated that equation (7?6o-) is of 
questionable value for pressure gradient flows.. Secondly, it has-been 
stated by the author that equation {'J.Sk) agrees well with experimental 
values of A* for non-zero y^. However, equations ,(7-6o) and (7*64) 
were both obtained from equation (’7*57) by utilizing the same set of 
assumptions. Consequently, the question arises as to how can 
equation (7*6o) be of questionable value, and. equation (7*^^) be of value? 
The answer lies not so much ia the derivation but more so in the response 
of the turbulent boundary layer to pressure gradients and changes in wall 
boundary conditions. It is well known that pressure- gradients affect 
the mean profiles mainly in the outer region, but equation {j,6o) Is 
applied only in the inner region; in particular, only in the viscous' 
sublayer. However, the mean profiles in the inner, region are particularly 
sensitive to wall boundary conditions which is precisely where 
equation is applied. 


A* = 26v* 


I? expp y*) 
p* J 
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Outer fiegion Model 

The eddy viscosity in the outer region is based upon the Clauser 
(ref. 32) model. The ratio of the eddy viscosity to the' molecular 
viscosity can be expressed as follows: 



(7.65) 


where 6* is the incompressible displacement thickness 'defined as 
^inc 

* • ’ 

6* = / ® (1 .. F)dy* ■ (7.66) 

^inc o 

' _ . 

The use of 6^ as the scaling parameter for the mixing length is 

discussed by Mai se' and McDonald (ref. 33 ) • KSee also, Morkovin,. ref.. -3^.') 

The value of in equation ( 7 * 65 ) is taken to be O.OI68 as reported 

in reference 35 * However, in order to account .'for the intermittent' 

character of the outer layer flow equation (7.63)' 'must be modified by, 

' . i 

an intermittency factor first obtained by Klebanoff (ref. 3 S)'> that is 



p*y* 


Kg-TF 


'-inc 


(7.67) 


where the transverse intermittency factor r(y) is defined as 


7 = 


1 - 


erf 


ii*. 


£ 

2 



(7.68) 
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Matching Procedure 

The criteria used to determine the boundary between the- inner and 
outer regions is the continuity of eddy viscosity. A sketch of a typical 
eddy viscosity distribution is presented in figure 5* 



Figure Matching procedure for two-layer model. 
The matching procedure may then be formally Written as follows: 


-) = — (1 - exp 

[ 1/1 , ji* 


r 


l/2 


y* 




i&u* 


A* 


Idy*- 


0 5 y* S-yJ 


/ ( 7 . 69 ) 


€\ 
[LjO 


p^pU* 

^ S* 7 

b* KinC 


y^ > y* 

■'m 


The location of the boundary separating the two region's^ y^ as determined 

V * ' » 

A 

from the continuity of equations ( 7 * 69 )j that is,,, where 




(7. 70)- 



55 


Eddy Conductivity 

In, Section 7-1-2 the eddy, conductivity was ionmlated in temiE of 
a static turbulent Prandtl number and "the eddy yi-scosity (see eq.s. (7 9-)j 

f 1 ‘ 

(7-10), and (7- 11))- The two-layer model eddy viscosity relation 
previously discussed suggests that there should be two models for the 
static turbulent Prandtl number^ cr^. However,, current experimental 
data for cf-t are inconclusive. The incompressible data which are 
available for the outer region of pipe flows (ref. 57) and boundary 
layers (ref. 38 ) indicate that has a' value between 0.7 to 0.9* 

These data indicate that as the wall is approached cr-t reaches a 
maximum value near the wall and then drops off rapidly to a value 
between 0.5 and O .7 at the wall. For the case of con^ressible flow 
there is very little data available. For a Mach number of 5 it appears 
that assumes a value very near unity in the outer region and 

decreases to values on the order of 0.7 to 0.9 at the wall (ref. 39)* 
However, Botta (ref. 4o) found that d,^ may achieve values as high as 
2 as the wall is approached. Bushnell and Becfewith (ref. 31 ) conclude 
that in the absence of better data for turbulent Prandtl numbers, the 
choice of the correct d^ distribution with y/S must depend upon 
agreement between the calculated and experimental profiles and wall 
fluxes. In the present analysis the turbulent Prandtl nuniber will be 
assumed constant with a value of 0.9 unless specifically stated otherwise 
however, any other value or functional relationship could just as easily 


be used. 



54 


T- 4 . 2 Transformed Plane 

Since the governing equations are solved in the transformed plane 
it is necessary to transform the eddy viscosity relations from the real 
plane to the transfomed plane. 

Inner Region 

In the inner region the ratio of eddy viscosity to molecular 
viscosity is as follows 


\ij^ i©5Sri 


(7.71) 


where Y is defined hy equation (7'. 44). The parameter JI appearing* 
in equation (7-7i) is the damping term and is defined as 


n = 1 - exp(- ]I-j_ Ilg) 


(7.72). 


where 


H - — 

^ \) T \© / 


{7.75) 


JI2 = 


^Pe^ j 0 b^e^o^w / 6 f ' 


26 I sJW V^. 


(7-74) 
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Outer Region 

■« ’ i 

In the outer region the ratio of eddy viscosity to molecular 
viscosity is as follows: ' ' ‘ . 




(T.75) 


where 


1 - erf 5(J/Yq - 0.78) 


(7.76) . 


and 
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'e 4--d 


t "'©(1 - F)dT) 


(7.77) 


7*5 The Transition' Region 

Equations (7-38).^ (7*39) ^ (7-^0) ;r (7-^2) together with the 
boundary conditions (eqs. (7*^5))r the' eddy viscosity relations 
defined by equations (7- 71) and' (7*75) complete the required system for 
either laminar or fully developed turbulent boundary layer flows. 
However^ the main, objective of the present analysis is to present a 
technique that will efficiently solve the laminar, transitional, or 
turbulent boundary layer equations as the boundary -layer develops along 
the surface. ‘Consequently, the location of transition the extent 

of the transitional flow - ^,i-> characteristics of the 
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mean flow structure in the transition region mst he taken into 
consideration. 

The author in developing this section on transition flow has 
had to purposely omit many 'excellent references because of the 
massive bulk of information available. A number of the cited references 
present extensive reviews on the subject and should be consulted for 
more information if required. 

7.5*1 Stability and Transition 

The stability of laminar shear flows and the eventual transition 
from the laminar state to the turbulent state of motion has probably 
received as much concentrated attention as any other problem area 
in fluid mechanics. The continuing interest in this area of fluid 
mechanics, is due to the fundamental importance of being able bo make 
a realistic prediction of where and if transition will occur on a given 
configuration in a 'given environment . The influence of many parameter s 
which affect transition has been studied and documentedj however,, after 
nearly sixty years of research contradictions or paradoxes still 
remain unanswered even for incompressible flows. Pioneers in transition 

♦ ’ i 

research were impressed by the three-dimensionality of turbulence and 
did not seriously consider the possible role of any two-dimensional 

^ 7 

amplification process such a^ predicted by the early stability work 
(see ref. 2). This attitude was imderstandable since the results from 
the- Tollmien-Schlichting (hereafter referred^ to as T-S) school of 
thought predicted an orderly vorticity pattern which differed 

1 

drastically from the observed random three-dimensional vorticity 
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components of transitional and ttnrbiolent boundary layer flows., 

- • m 

^ _ 

problem was further complicated by the high-level of free- stream 
■turbulence that existed in the early wind timnels which apparently’ 

*■ * t 

prevented the development of the’ T-S -type amplification process. 
Furthermore j due to the limited amount of transition' data and its 
apparent agreement with Taylor's (ref. 4-1) local .separation- criteria 

f 

there was no widespread acceptance of the T-S process xmtil about 1947* 

i 

•' f ' 

The first indication that the T~S amplification process existed was 
obtained by Schubauer and Skramstad (ref. 42.) in their classical 
-experimental study of flat plate flow in the low free-stream 
turbulence tunnel of the National Bureau of Standards. 

The T-S instability arises- from the. viscous instability in the 
laminar shear layer. It is a two-dimensional wave motion that 
selectively amplifies disturbances over a band of unstable frequencies. 
C. C. Lin (refs. 43 and 44) was able to calculate with very good 
agreement the results of the- Schubauer-Skramstad experiment . Klebanoff 
and Tidstrom (ref. 45) were able to trace the growth of the unstable 
. waves into the nonlinear region -of transition where the T-S approach 
becomes invalid. They observed that when the velocity perturbation 
due to the wave exceeded approximately 10 percent of the meein .speed the 
nonlinear effects became important. Any further development of the 
mean flow instability was- then three-dimensional. In this region of the 
transition process "peaks and valleys" form consisting of variations in 
mean cross flow. The "peaks and valleys" eventually break down and form 
turbulent "spots." 
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A turbulent spot is essentially an isolated region of turbulent 
flow surrounded an otherwise laminar flow: field. It, has. been estab- 
lished that the flow within the spdtihas'the same basic. characteristics- ’ 
as that, of a fully developed turbulent flow. ‘The existence of these 
isolated spots of turbulence was first reported, by Emmons (ref. 46) 
in 1951 (see also ref. 47)* KLebanoff, Schubauer, and co-workers 
(ref. 48) were the first to carefully study the shape of the spots and 
■their characteristics. They confirmed that the flow within a given 
spot was turbulent in nature , that the spot propagated through the 
surrounding flow field., and that the spot grew in size as it was 

1. 

convected downstream. The transition process is. completed by either, 
or both, the transverse growth of spots or cross-contamination and the 
creation or production of new spots at new transverse locations. When 
the entire flow region at some streamwlse location is covered 

ly a turbulent flow structure the boundary layer is then said to be 
fully turbulent. 

The natural transition process ly which a fully developed flow 
structure is obtained may be categorically di"vided into six separate 
stages. The first stage, of course, is the laminar region of the 
boundary layer. The second stage begins with the formation of two- 
dimensional unstable waves arising from an instability in the laminar 
layer. The third stage, which marks the beginning of the nonlinearity, 
is the formation of "peaks and valleys" or longitudinal streaks in the 
mean velocity profiles. The fourth stage occurs with ■fche breakdown 
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of one or more of the low- speed streaks and the. formation of "-spots" 

i 

of turbulence. The fifth stage consists of the growth of -the spots 
through cross contamination and the process of .new spot creation 
(production) • The sixth and final stage in natural transition is then 
fully developed turbulent layer. (See Kline j S.J., reference 49.)' 

Stability theory cannot currently be 'used to predict either the 
nonlinear details of the transition process after the two-dimensional 
waves have been -amplified or the location of transition, x^. j_. 

Stability theory can,, however, establish which boundary layer profiles 
are unstable and the initial amplification rates. The theory can 
identify those frequencies which will be amplified at the greatest rate 
as well as. the effect on stability of various flow parameters. One of 
the more important contributions of linear stability theory has been 
to reveal -how these flow parameters should be varied in ozder tO' delay 
transition, i.e., wall cooling, suction, etc. The Blasius profile was 
studied by Tollmein (ref. 50 ) in 1929* The results of his work remained 
xmconfirmed experimentally until the classical experiments of Schubauer 
and Skramstad in 194? (ref. 42). Since the beginning of the past decade 
the .solution techniques utilized in stability theory have been greatly 
modified by the availability of high-speed digital computer systems. 

Now, instead of having to work many hours in order to obtain a minimum 
-of results of questionable accuracy, the digital ■ computer can obtain 
an immense quantity of highly accurate numerical resul-fcs from the 
governing system of equations. 'A re-yiew of methods used to predict the 
location of transition, from stability theory is- presented by 'jaffe. 
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Okamuraj and Smith, in reference 51* However, a^poi'nt which should he 
strongly stressed is that a thorough study of the connection between 
stability and transition still remains to be completed. 

There currently exists a number of good review articles and 
books on stability theory. Review articles have been written by Stuart 
(ref. il-), Shen (ref. 52), Reed (ref. 53), and Drazin and Howard (ref. 5^)*- 
The article by Reed considers asymptotic theory in detail. Drazin and 
Howard in turn consider only the inviscid theory. The articles by 
Schlichting (ref. 2) should also not be overlooked. The most complete 
study of modem boundary- layer, stability, in the opinion of the author, 
is presented by Mack (ref. 55). , 

7 . 5*2 Transition Location 

Many parameters influence the location of transition. These 
parameters can best.be thought of as forming a parameter phase space. 

Such a parameter phase space would include Reynolds number, Mach number, 
unit Reynolds number, surface roughness, nose bluntness, pressure 
gradients, boundary conditions at the wall, angle -of -attack, free- . 
stream turbulence level, and radiated aerodynamic noise. Morkovin 
(ref. 9 ) recently con^jleted a very extensive and thorough examination 
of the current state-of-the-art of transition in shear layers from 
laminar to turbulent flow. The most striking conclusion that one 
obtains from the review is that although d great hulk of experimental 

/ K • » 

data on transition currently exists, much of the information on high- 
speed transition has not been 'documented in s'ufficient detail to allow 
the separation of the effects of multiple parameters on transition. A 



discussion of the effects of eacji of the parameters that may influence 

I * , 

transition in- high-speed flow is beyond, the scope of the present paper. 
Furthermore, It is seldom possible to study any one parameter experi- 
mentally while holding the remaining parameter phase space constant. , 

The reader, interested in a -detailed discussion^ is directed tp the - 
paper by Morkovin (ref. 9') where over 500 related references are cited 
and discussed. The paper by Savulescu (ref. 8); should also’ t)e consulted. 
Another,, although less' detailed discussion, is- presented by Fischer 
(ref. 56). The effects of radiated aerodynamic noise on transition is 
discussed by Pate and Schueler (ref. 57) • -Hypersonic transition, to 
name but a few references, is discussed by Softley, Grabel, and Zemple 
(ref.. 58),- Richards (ref. ^ 9 ), Potter and 'Whitfield (ref. 60), and 
Deem and Murphy' (ref. 61) . A very good discussion on the effects of 
extreme surface cooling on hypersonic flat-plate transition is presented 
by Cary (ref. 62). The effects of nose bluntness and surface roughness 
on boundary layer transition are discussed by Potter and Whitfield 
(ref. 63). 

It would be prohibitively time consuming to obtain a detailed 
stability analysis for each of the test cases considered in the present 
paper (see Chapter DC). Furthermore, even if such an analysis were 
completed the transition location would not be obtained. Consequently,, 
in the present analysis the location of transition will be determined 
by one of, or a combination of, the following three methods. These 
methods are (l) a stability index or vorticity Reynolds number first 
proposed by Rouse (ref. 64), (2) correlations based upon a collection 
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of experimental data over a 'broad range of test conditions', ^ and ( 5 ) i>y 
using the meastired experimental id’cation of transition as a, direct 
input into the analytical solution. One might immediately a'ssume that 
the latter of the three methods would be the most accurate j however, 
caution must be exercised in interpreting and applying the results, of 
experimental transition measurements since data on all parameters that 
may affect transition are seldom available . Also-, the transition 
location may very well he a function of the method used to obtain the 
•experimental data. There are a number of ways to experimentally 
measure the transition location which often yields different locations 
for the identical boundary layer. This fact has contributed to much of 
the current confusion and scatter that results when attempts are made 
to correlate a great number of experimental transition data obtained 
in a number of different investigations without considering how transition 
was determined. 

7 . 5*2.1 Stability Index 

Hunter Rouse (ref. 64) nearly 25 years ago obtained through a 
process of dimensional analysis combined with observed fluid motion that 
.a stability index could be expressed as follows: 


Su* 

V* By-Sf- 


(7.78) 


This index has the form of a vosrbiclty Reynolds number which is 
obtained from the ratio of the local inertial stress •p'^*^( 8 u*/Sy*)^, 
to the local viscous stress ii-5f(Su*/^*) . Rouse assumed that in order 



for transition to occur the stability index should reach some limiting 
value which was assumed invariant. He was able to further show that 
this invariant value, ' (^max) cr should be on the order of 5D0 for 
incompressible flows. . , ’ * 

The use of ^ stability index is- in principle' similar 

to the basic reasoning which led'O. Reynolds (ref.' 7) in I885 to 
postulate that the nondimen sional parameter ud'/v • could 'be used to 
define a critical value (ud/v)^j, at which transition would occur in 
a circular pipe of diameter, d. 'Unfortunately, (Xjji£Q5;)cr function 

of the transition parameter phase space in much the same fashion as the 
critical Reynolds number, and cannot in reality be a true invariant 
of the flow as suggested by Rouse ;(see Section 9-8j fig. 19). The 
stability index does, however, possess a number of important character- 
istics which can be .directly related to the stability of laminar flowsj 
for example, the position of the critical layer can be obtained directly 
from equation (7.78). 

A typical transverse distribution of X for a compressible 
laminar boundary layer is presented in figure 6. The stability index 
has a zero value at the wall and approaches zero as the outer edge of 
the layer is approached. The maximum value of X, X mg-v > will occur 

(y/S) are of importance in the usage of the stability index as a 

^ax 

guide to boundary layer stability and transition.' 


at some transverse location, (y,/B) 


^ax 


The values of X 



Ifif 



Figure 6.- Vorticity Reynolds nuinber. 

As the laminar boundary layer develops over a surface, X 

max 

increases monotonically until the crftical value (^max^cr reached 

at which point transition is assumed to occurs that is, the location of 

X .. For compressible flows (x ) is "not an invariant. In the 

’'Cz* 

present study was found to vary frcm approximately 2100 to 

values on the order of ^1-000. The variation of is a strong 

ftinction of unit Reynolds number for data obtained in air wind tunnel 
facilities (see Section 9-®) s-s would be expected from the previously 
mentioned transition studies. However, while not invariant the 
stability index does exhibit the same dependence on various parameters 
as results from the more complicated stability theory. For example, at 



is found' to decrease 


a given streamwise location, *, x,. the value.. of 
(which implies a more stable flow), with wall cooling, wall suction, and 
favora'ble pressure gradients, whereas it increases (.which implies ,a more 
unstable flow) with wall heating, mass injection ("transpiration), and 

’ , ' i 

adverse pressure gradients. To the author’s knowledge the stability 
index has been used as a correlation parameter in only two boundary- 
layer transition studies. A modified form of the -parameter was used to 
correlate the effect of free-stream turbulence on transition by Van Driest 
and Blumer (ref. •65)* Correlation attempts, using Rouse's original 
invariant assumptions, were made in reference 66j however, the results 
were only fair. 

One of the most important characteristics of the vorticity 

Reynolds number is that the value of (y/S)y is in excellent 

'^max 

agreement with the experimental location of the critical layer which 

represents the. distance normal to the wall at which the first high 

frequency burst representing laminar flow breakdown will occur. 

■Stainback (ref. 67) recently Obtained the Rouse stability index for 

similar laminar bo-undary layer flows over a broad range of wall-to-total 

temperature ratios for Mach numbers up to 16. The numerical calculations 

were made for both air and helium boundary layers. The agreement between 

(y/S)y and ‘the experimental critieai layer position was excellent 
'hnax 

over the entire range. ('See Section ) 
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7-. 5 • 2 . 2 Empirical Correlations 

In most instances the design engineer has to rely on empirical 

correlations of experimental transition data in order to fix the most 

prohahle transition location for a design configuration. However, 

caution should always he used when obtaining the most prohahle location 

of transition from such correlations, since' any given correlation is 

based upon a specific 'collection of data which will not he completely 

general. Furthermore, the transition process may not be unique;- for 

example, the blunt-hody transition paradox (see refs. 9 and' 68). 

* * « 

There currently exists a large number of empirical .correlations 

1 .* 

for predicting the p.rohahle location' of transition. Some of these 

* > 

correlations are of questionable value; however, some can be used with 
confidence providing it is realized that one is predicting' a probable 
range of locations- and not an exact fixed point. - One of the more 
successful correlations was' obtained, by Beckwith .(reft - 69 ) -at the 
Langley Eesearch Center. The correlation developed by Beckwith is 
based on experimental transition data obtained over a wide range of 
test conditions in air wind tunnels, ballistic ranges, and free flight. 
The correlation can be expressed as follows:' 

0.7exp(-0.05Mf) 

= I + JM (0 ) 
e w 

where E* denotes the unit Reynolds number per foot, vl*/v*. The 
constants I and J are functions the environment in which transition 
was measured and are given in the following table. 


{7.79) 


log 


10 


Re 


St 
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Facility 

I 

J 

Air wind tunnel 

0.95 

0 . 167 ’ 

Ballistic range 

1.00 

0.125“ 

Free flight 

1.52 

0.150 , 


Equation (7-79) .can be expressed in terms of the transition Reynolds 
number, follows (see ref. 69 ): 


Re- 


fa 

10 L 


J^,*0.6 _|2I + 2JMe(©w) 


2 — I 

0.7exp(“0.05Mg) 


2 I 


w 


(0.094 M| + 1.22 0^)^ 


(7.80) 


For an in-deptb review of the problems associated with transition the 

> 

reader is referred to the results of the Boundary Layer Transition Study 

* t - < 

Group meeting held in San Bernardino-, California,, in 196'J'»(refs. 69 , . ' 

70, 71; and- 72) . . 

7 . 5-2. 5 Experimental Transition 

Much of the confusion that exists today concerning boundary 
layer transition may be attributed to -one of, or a' combination of, the ' 

' t 

* - L ' ' 

following three factors.- The first factor is that in many instances the 
investigator who made the experimental stu^y may not 'have carefully 
measured or recorded the exact conditions under which 'the data were • 
obtained. The second factor is that the experimentally observed transitior 
location depends on the experimental technique used to obtain its 
location. The third, the most pessimistic factor, is that transition 
may not be a unique process under all conditions (ref. 9) • 
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The importance of carefully measuring the environment under 
which the experiments are made cannot he overstressed'. In the past, 
many of the factors 'which may influence .transition such as free- stream 
turbulence and- acoustic radiation from the tunnel side' wall boundary • 

' I ' . 

layer were not measured. (See references 69 to' 72.) The location of 
transition as obtained experimentally is a strong function -of the method 
used to- determine its location. There, are currently a number of 
techniques used to obtain the transition location. , Some of these 
methods are hot-'wire .traverses, pitot- tube surveys near' the wall, 

1 '' 

visual indication from- schlieren photographs', and heat transfer 

measurements at the wall,.^ -Each of these methods basically measures a 

different flow process. Consequently,, it- would be misleading to believe 

* ♦ 

that each technique would yield the same location for transition if 
simultaneously applied to the same boundary layer. Of course, the 
concept of a transition '"point” is misleading in itself since transition 
does not occur at a "point" but instead-over some finite size region. 

For the test cases .presented in the present analysis the 
experimental transition location will -be determined from heat transfer 
measurements at the wall whenever possible. The main reason for this 
choice is that it is the method most often used in the literature. 
However, .it should be noted that the actual nonlinear transition 
process begins somewhat upstream of the location where the heat transfer 


at the wall deviates from the laminar trend. 



7.5*3 Transitional Flow Structure 


Assumng that the transition location has been fixed for a 
given problem one must next consider the following two important 
factors; firsts the length of the transition regiouj 
(sometimes referred to as transition extent)^ and secondly^ the mean 
flow characteristics within the region. Once appropriate models, are 
obtained for- these two factors it is possible to smoothly connect all 
three flow regions such that one set of goveming equations may be 
used. 

The classical laminar boundary layer equations should' yield 
reasonably accurate profiles and wall fluxes in the linear region of 
transition; that is^ prior to the turbulent spot formation. The 
intermittent appearance of the turbulent spots and the process of 
cross- contamination is not well under stoo*d..' The spots originate in a 

more or less random fashion and merge 'with one another, as they grow 

.. * 

I 

and move downstieam. Eventually, -the entire layer is contaminated 
which marks the end of the transition process^ process by 

which cross-contamination occurs appears to haye been' studied in detail 

' f 

only by Schubauer and Klebanoff (ref. 75) since its d'iscoye'ry by 
Charters (ref. 7^) la 19^5 • As the turbulent spots move over a fixed 
point in the transition region, the point experience s ' ah alternation 
of fully laminar flow when ,no spot is present to fully turbulent flow 
when engulfed by a spot. These alternations can be described by an 
intermittency factor which represents the fraction of time that any 
point in the transition region is engulfed by fully turbulent flow. 
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The distribution of spots in time and space in Gaussian for , 
low speedy natural transition. However^ very little is known about, 
the spot distribution in hi^ speed compressible flow. Furthermore^ 
there is no assurance that the spot formation and distribution in 
hypersonic flows will be analogous to the low speed model . -However , 
in the absence of a more satisfactory theory, the -author has chosen" 
the approach of Dhwan and Narasimah (ref. 75) which was obtained 
mainly from low speed data. In reference 75 the source density 
function of Emmonds (ref. k6) was used to obtain the probability 
distribution ( intermittency ) of the turbulent spots. Dhwan and 
Narasimah (ref. 75) obtained an intermittency distribution of the form 

P(I) = 1 - exp(- 0.412 p) (T-8l) 


where 



(7.82) 


for < X ^ ^,f • term | in equation (7«8l) represents a 

normalized streamwise coordinate in the transition zone, and A is a 
measure of the extent of the transition region) that is 

A = (x) - (x)* ^ ' (7-85) 

r'=3/4 • r=i/4 

' S. ' t 

In the present analysis, the tr^sition e^ent, Xj. ^ ^1^ (see 

‘ . ‘ ^ . 

Section 7*5.4) is first calculated; next, the intemittency function 

* \ ^ ^ r , 

• i 

is determined such that P = 0-9-999' at x = " 
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For low speed transition, 'Dhawan and Narasimah show that all 
transition regions may he correlated to form a universal intermittency 
distribution. Potter and Whitfield (ref. 65 ) have shown that similarity 
exists in the transition process for subsonic, supersonic and hypersonic 
flows. Althovigh the effect of pressure gradients on the intermittency 
distribution has not been studied in detail, it is believed that the 
distribution will be affected primarily in the region of the initial 
spot generation. 

The intermittency distribution- across the layer at a given 
streamwise location in the transition region is a function of the 
shape of the turbulent spots. In reference 75 it is shown that the 
spots have a nearly- constant cross sectional area close to the surface. 
The intermittency distributions in the transverse direction (y-direction) 
are similar to those observed by Corrsin and Kistler (ref. y6) for fully 
developed turbulent boundary -layers (see ref. 77)- Corrsin and Kistler 
found that the transverse intermittency varied from a maximum of unity 
near the wall to a near zero value at the outer edge of the layer. The 
transverse intermittency distribution is of secondary importance in 
relation to the streamwise distribution ‘in determining the mean profiles 
and wa-ll fluxes. In the present an^ysis the only intermittency 
distribution applied in the transverse*, direction (y-direction) is that 
of Klebanoff (ref. ^6)- in the outer layer as applied to the fully 
developed t-urbulent layerj that is, equation (j-.:68'). 
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7-5*^ Transition Extent 

The assumption of a universal intermittency distribution in^jlies 
that the transition zone length (transition extent) can he expressed 
as a function of the transition Reynolds number, i-/'^e’ 

reference 75 it is shown, for the transition data considered, that the 
data are represented on the average by the equation 


= ctRe^ (7.8^) 

■^,i 

where Ax^ = The coefficients a and p are found to 

assume -values of 5-0 and 0.8, respectively. The location of the 

completion of transition, 'xf' „ can then be obtained directly from 

■ 6,1 

equation (7.8ii-) as follows: 

+ oRr^(Re^ y (7.85) 

where Re-*- is the local unit Reynolds number, vl*/v*. . Morkovin 
(ref. 9') found that only about 50 percent of the experimental data he 

j- 

•> 

considered could be fitted to the low speed universal curve- of Dhawan 

V * 

and .Narasimah; that is to equation (7.8^+). This was .to be expected, 
since the data considered in reference 75 covered only a very limited 
Mach number' range . ' . - * ^ , 

Potter and Whitfield (ref. 63) measured -the -extent of the 

i ‘ ^ 

transition zone over a rather broad Mach number' range (3 — 5.; 

1“^ =8).' They observed that the transition region^ -when defined in 
terms of Re^.^. is basically independent' of the -unit Reynolds number’ 
and leading edge geometry; that is 
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( 7 . 86 ) 


They noted (ref. 65) that the extent of the transit;ion region increased 
with increasing transition Reynolds nimiber over the Mach member range 
0 < < 8 for adiabatic walls. The extent of the transition 

region was also observed to increase with increasing Mach numbers for 
a fixed transition Reynolds number. 

In the present analysis^ due to the. lack of general correlations 
for the extent of transition this quantity will be obtained directly 
from the experimental data unless otherwise noted. In particular, if 
heat transfer data are available the transition zone will be assumed 
to lie between the initial deviation from the laminar heat transfer 
distribution and the final peak heating location. The transition 
region defined on the basis ,of the Stanton number distribution is 
presented in figure 7 . design engineer does not have the advantage 

of experimental data which were obtained under the actual flight 
conditions. Consequently, the most probable location of transition 
would be obtained from a correlation such as presented in equation ( 7 . 79 ) 
and (7.80). The extent of transition could then be obtained from 
equation (7*85) or an approximate relation such as follows; 

c 2 (7.87) 

■ 

Since the main objective of the design engineer is to ob.tain esti^tes 
of the heat, protection requirements or viscous drag, for example, the 

i ' 

errors involved in using equation (7.87) would be negligible for a ^full 



Stanton number, Ch 





Figure J,- Transition extent definition. 
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scale vehicle. However^ inaccurate estimates of >the transition 

' ' T 

location^ ^ could lead to significant errors. i 

^ 

j * 

* f 

7.6 Boundary Layer Ea-rameters 

Certain lioundary layer parameters are of interest to the design 
engineer as well as to those interested in the general development of 
the houndary' layer for arbitrary boundary .conditions. A number of 
these parameters', which are used in chapter IX, are presented in this 
section for the physical plane as ■well as for the transformed plane 
where meaningful . 

7 - 6.1 Displac^ent Thiclmess 

The displacement thickness is a measure of the distance by 
which the external streamlines are shifted because of the' development 
of the viscous boundary layer. If in the problem under consideration 
interaction is in^ortant, .the displacement thickness may be used to 
determine a new effective body. (ref. 15) . The boundary layer dispilacement 
thickness for compressible flows is defined as 

a* = tj(i - (7.88) 

Jq Pe"^ 

where 5^* is defined as the nondimensional displacement thickness, • 

S*/l*. In the transformed plane the nondimensional displacement 
thic3mess becomes 

-H* = — - F)dri 


(7.89) 
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It should be noted that the definition of 5* for axi symmetric 
flows with transverse curvature is not the same as the definition of the 
nondimen sional displacement thickness for two-dimensional flows^ ^*2D’ 
The displacement thickness with transverse curvature is related to its 
two-dimensional counterpart by the following relation (ref. 78 ) 


r 

o 


26, 




2D 


( 7 . 90 ) 


The incompressible displacement .thickness used in the outer 
region model for the eddy viscosity is presented' in equation (7*66) and 

4 

( 7 . 77 ) for the real and transfomed pr^esj "respectively. 

7 . 6.2 Momentum Thickness ' r - « 

' J- ^ L ^ ' 

The momentum thickness is. tised in a number of solution techniques 

. 1 . ■ 'i ' 

. I ' ' ♦ - 

as well as. in current transition correlations.' The' nondimen sional 


momentum thickness for the physical plane is defined as follows; 


Jq Pe^e V 


( 7 . 91 ) 


The nondimen sional momentum, thickness in the transformed plane may 
then be' expressed as follows ; 


0 


F(1 - F)dll 

Pe^e^o ° 


(7.92) 



T.6-5 Shear Stress 


The viscous drag of a body moving through a fluid is dependent 
upon the shear stress at the wall. For turbulent boundary layers the 
shear stress distribution across the layer is also an important 
parameter. The total nondimensional shear stress may be' divided into 
two components, laminar and turbulent; that is 


^t 


(7-95) 


or,, in terms of the current notation 


- 5u 

= pe — 

oy 

•The nondimensional total shear stress. Tip is defined as 


( 7 . 9 ^) 


TT 'fp/ 


r^ ref^ef\ 

L* } 


(7.95). 


Equation ('7*95) becomes in the transfojmied plane 


2 ij. - 
Pehe^e^o"^ 1 e 


Trp - 


os 


ViF n, 


( 7 . 96 ) 


where e is defined in equation ( 7 - 15 ) • 

As the wall is approached, Tp approaches and equation (7-96) 

becomes at the wall 


PgPgUgrQ / 

Tw = —{^ — 


( 7 . 97 ) 
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7‘6.k Skin Friction Coefflcxent 

The skin friction coefficient is usually of more interest to the 
design engineer than the actual magnitude of the shear stress. The 
friction coefficient is defined as 


- ■* / 1 * 
Cfg = ■^T/ 1,2 


(T-98) 


which becomes in the nondimensional physical plane 


o 2' 

■2oi T, 


T 


2 

Pe^e 


(7.99) 


or 


Cf. 




Pe^e 


2\^yJ' 


(7- 100) 


In the transformed plane equation (7-100) become s* 




(7-101) 


or^ when evaluated at the wall 


2qie^o / 


I r- 


\fir \ 


w 


(7.102) 
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T»6.5 Heat Transfer 

The heat transfer between the boundary layer and the wall is 
required for aerodynamic heating analysis. The heat transfer at the 
wall is related to the temperature gradient at the wall in the 
dimensional physical plane through the following equation: 

(7.105) 


or in nondimen sional form' as 


^Lw = - - 


ii ■ 

Ww 


(7.1o4) 


where 






(7.105) 


Equation (7-10^) becomes in the transformed plane 


•Iw = - 


W^eVoA 


( 7 . 106 ) 


7 . 6.6 Heat Transfer Coefficient 

The heat transfer coefficient is defined by the' following 
relation: 


h* = c||/(TJ - T*,) 


(7. '107)' 



6o 


where T* denotes the dimensional adiabatic wall temperature. In 
aw 

the present analysis the adiabatic wall temperature is defined as 


^Iw = (T.IO8) 

Equation (7*107) can be written in nondimensional form as follows: 


h = 


^w ” ®aw 


(7.109) 


where the nondimensional heat transfer coefficient h is defined as 


h = h* 




L* 


'■ (7.110), 


The recovery factor r^ used in equation '(7.108) has" the .following 


form 


Tf =\/~^[l + (cr“^/^ 1)1^. 


(7.111) 


This relation was developed by the author and 'is. based' dn the;concept 

* » ^ . * 

that the flow within the turbulent spot has ’the -same structure as that 
of the fully turbulent region of the boundary layer. ,It can be -seen 
from equation (7. Ill)' that, the recovery factor assumes the accepted 

1/2 

mean values for pure laminar or fully turbulent flows of .o'' and 

respectively. Within the transition region the recovery factor 
is a function of the turbulent spot distribution through the inter- 
mittency factor,, E. (See equation (7.81).) 
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7 . 6.7 Stanton .Musiber 

A convenient nondimensional parameter often used in the 
literature to. describe the heat transfer at the -wall is the Stanton 
number which is defined as follows: 


Cl, =h*/(cjp*u*) (7.112) 

e 

7*6.8 Friction Velocity 

The velocity scale most often used for turbulent flows is the 
friction velocity. This important parameter is defined as' 



which becomes in the transformed plane 


(7.115) 


u, 


T- A 


r 5f\ 
r ^ 

\ /w . 


\iw 


(7.1lit-.> 


where 


' ' . * /V* 

,u = ,u /u 




T T ref 

7 . 6.9 Law of Wall Coordinates . ” , • , 

The iaw-of~the--wall coordinates used in 'turbulent flow 
(see eqs. (7*^) and (7*^9)) ane defined' 'in the" physical pl^e as 


(7-. 115)' 


y = y 


* * / * 
V 


(7.116a) 
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and 


u =. u 




(r.ll6b;) 


Equations (7-H6-) may be written in the , transformed variables as 


• = 


Pe^T~^ 

cc(i 2 
© 


■(7.117a) 


and 


u+ = F- 




0 


- 1 - 1/2 


V 


(7.117b) 


7.6.10 Defect Law Coordinates 

The defect law coordinates used in turbulent flow -(see 
eq. (7.50)’) are y/b and u. The coordinate u is defined as 
follows; 


u = (u 


* 



(7.118) 


or 



(7.119) 


7.7 Fitot Pressure Profiles 

The total pressure is one of the mean profile values that can 
be measured directly in experimental studies. -Consequently, it is 
advantageous to calculate the total pressure distribution across the 
layer for -comparison to experimental data when available. The total 
pressure distribution is often encountered in the experimental 
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literature in one of tvo forms; (l) the local isentropic total pressure 
distribution referenced to the isentropic value of the total pressure 
at the boundary-layer edge (the term isentropic as used here means that 
the flow is locally brought to rest isentrppically)', and ,(2) the actual 
pitot tube measurement of the total pressure referenced to the pitot 
tube value at the edge of the bounaary layer. For subsonic flows the 
two methods are identical; however, for supersonic flows caution must 
be exercised in interpreting the experimental results. 

If the static pressure is constant across the boundary layer, 
the local isentropic total pressure referenced to the isentropic value 
at the edge of the layer may be expressed as 

r/r-i 


Pt _ f2 + (7 - 
Ptg L2 + (7 - 


( 7 . 120 ) 


The local Mach nmber, M, may be expressed in terms of the mean 
temperature and velocity profiles as followh': 


M = i . ^ (7.121) 

\/(7,- i)T^e 


which becomes- at the outer edge . ^ 

‘ * A 

(7.122) 

\/(7 - l)Te 

' * * 

In the case of subsonic flows the experimental profiles can be 
directly compared to the distribution obtained from equation (7.120); 



however, for supersonic flows it is necessary to divide the boundary 

t 

^ f 

layer into a sub son i c and a supersonic region. '[Uhe boun^ry separating, 
these two regions occurs where- the local Mach nvariber'is unity and is 
obtained from equation (7*121) i , In the subsonic region of the boundary 
layer the pitot tube measures the -local isentropic total pressure; 

however, it is usually referenced to the measured pitot value at the 

‘ ■ ) ' , - ' 

. > .* 

•edge of the boundary layer where the flow is supersonic. ' Then the 
correct expression for the sub-sonic region- becomes ■ ' ‘ 






'2 + (r -r i)m^ 


(7 + i ) m 2 




7 + 1 


1/7-1 


(7.125) 


In the supersonic region of the layer the ratio can be written as , 


27/7-1 


r 


“m 






27M^ - (7-1) ' 

27M| - (7 - 1)' 


-,-1/7-1. 


(7.124) 


In the present analysis the total pressure distributions will be 
obtained from eq.uations (7. 123) and (7.124) for the subsonic region and 
supersonic region of the boundary layer, respectively, The boundary 
between the two regions will be obtained from equation (7.121); that is, 
where the local Mach 'number is unity. 


7.8 Exterior Flow 

The following quantities will be specified in order to describe 


the exterior flow and fluid properties 
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Reference Values 

TRe reference values are evaluated for an isentropic^ perfect 
gas as follows: 



( 7 . 125 ) 


2ref 




7-1 


.-1/7-1 


(7 - D'cJt* , 


.(7.126) 


(T-I2T) 

2.270(T* )5/2 X 10-*-^‘ 

M. J ; (air only) , ( 7 -^ 28 ) 

^^?ef + S^ ■/ ' . • 

The reference values are independent , of hodjr geometry and' are used to 
form the dimensionless variables (see eq.s. (7.'20’)).' 

Edge Values ’ - 

^ * ** V J V “ * 

4 * * 

■ The inviscid pressure distri'bution must "be specified , for a 

I \ * 

given flow geometry in qrder to obtain a solution of the governing 
differential equations (eq_s ./ ( 7-..38) to, (7;lo)). This distribution can 

‘ r 

t. be obtained from either experimental data or from an exact solution of 
the full inviscid Euler equations. In the case of supersonic flows 
where shock wave curvature may have a first order effect on the 
boundary layer structure^ the effect of variable entropy should be 
taken into consideration (see Section 7*9) • The total pressure behind 
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the shock wave , can he obtained from' the oblique shock relations 

(ref. 79) provided the shock wave, angle^ 9 g(rg,Zg) is specified; 
that is j . • • ^ 



(7+ 1)M^ sin^0 (r ,z ) 



> 7-1 


7 + 1 ' ■ 


(7 - sin20g(rg,Zg) + (271^ sin 26 g(rg,Zg) - (7 - 1 ) 


(T-129) 


where is the dimensionless total pres stir e behind the shock wave 

s 

evaluated locally at the shock wave coordinate (r ,z ). For 'cases 

3 5 

where the effect of variable entropy is either neglected or -negligible-, 

P^ is evaluated at r^ = 0 (see fig. 1) and is invariant along the 
s 

streamline i|fg(o,x). For cases where the flow is either subsonic 
everywhere or supersonic shock free flow, P. = P4. 

■Cg -c,co 

IPhe edge conditions can then be calculated from the specified 

static pressure distribution and the known value of Pt^ as follows 

s 


Mg = 



7-1 

7 



(7.130) 



(7.131) 


u„ 


= MgV7; 


r)T. 


(7.152) 


p 


e 


7 



(T-133) 



(7-13^) 


^a.= (Ta)^^^|^). (airily) 

^ ' 

where S = . ' ; - ■ ‘ 

t ^ ' , 

V 

7*9 Va-rlable Entropy 
One of the important factors that should be taken into con- 
sideration for slender blunted or sharp-tipped bodies in supersonic, 
flow is the effect of the variable entropy introduced by shock wave . 
curvature on the boundary layer characteristics, ihe basic problems 
associated with incorporating variable entropy into the boundary 
layer solution are (l)’ to obtain the entropy distribution alo33g the' 
boundary layer edge as a function of and (2.) to develop an 
iterative solution technique with which this distribution can be 
efficiently included in the solution of the governing boundary layer 
equations (eqs. ( 7 . 38 ) to (7.4o)). In the present analysis the shock 
wave coordinate's and inviscid pressure distribution are obtained from 
an exact inviscid flow field solution. The boundary layer equations 
are then solved by the method presented in chapter VUI. At each x- 
solution station a mass balance is obtained between the total mass 
contained in the boimdary layer at the particular x-station and the • 
mass passing through a streamtube of radius r* (see fig. l) . Ih'om ' 
this the entropy distribution as a function of x can be obtained. 
The resulting, total pressure at each x-staticn is stored for future 
use in the digital computer. The actual equations and procedures 
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utilized in the process a3?e discussed in the remaining porticn of this 
section . 


The shock wave coordinates and inviscid pressure distribution 
can he directly obtained from the solution developed hy lomax and Inouye 
(ref.. 8o) for -either, sharp-tipped or hlunt axi symmetric or planar 
geometries. The mass balance equation which equates the mass of the 
fluid passing through a streamtube of radius r* (see fig. l) to the 
total mass' entrained in the boundary layer at a given x- station, can be 
expressed for either axi symmetric or planar flows a-s follows 


p* u* ~ (a + 1) J dysf- 

ref ref s o 


(7-155) 


or in terms of the dimensionless variables as 


=,(-j H- ® I'rg dy •. 


(7-136) 


Equation (7-136) can. then -be expressed in the transformed plane as 
follows:' 


rg = 


I 

(j t lh\[H F 


dT)- 


jtl 


(7-137) 


Equation (7.137) is locally solved for rs at each x-station 
along the body surface. The shock wave angle,, C^^s-j^s) iken 
obtained from the known shock wave coordinates at the shock wave 
streamline intersection point, total pressure behind 
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the shock wave oa the streamline then evaluated from 

equation (7*129) • These values are stored' In the computer for each 
x-station. This process is continued until the solution at the final 
X- station (maximum x- value) is qompleted. The- houndaiy layer solution 
is then known for the_ condition of constant entropy. This is called 
the first iteration and if variable entropy is to be neglected it 
represents the final solution. In .order to include variable entropy 
effects the solution is. then restarted at | = G (x = 6) and the 
solution recalculated 'using the previously stored values of J*.}- 

* I » ' ’ s 

' , * • ‘ • * 

to -calculate the new local values 'of Mg^ Tg^ Ug^ Pg'^ and pg (see 
eqs.. (7*150') to (7*13^)).i th'isMs defined -as the second ite'ration. It 

i 

should be noted that Ej. is Independent of, | for the first 
iteration cycle; however, for subsequent iteration cycles is a 

function of g* The ne-^r values for the i + 1 iteration .are 

calculated .during the i iteration. The Pg distribution as a function 
of X is an invariant for -all iteration cycles. 

The convergence criteria used at each x-station is as follows: 




< Q 


(7.138) 


where i is the iteration index and Q denotes the specified 
convergence requirement. For engineering calculations three iterations 
are generally sufficient for a.- one percent or smaller change in Ug at 
each x-station; that is, for Q = 0.01. 



VIII. NUMERICAL SOLUTION QF THE GOVEKKTIWG- EQUATIONS 


The go-veming equations for the compressihle latninar, transitional, 
and turbulent boundary layers form a fifth-order system. Three of the 
equations are nonlinear partial differential equations (see eqs. (7*38-) 
to 7.Ji-0)) and the remaining two are algebraic relations (see eqs. (7*31) 
and (7.52)). The most important feature of this system is that it is 
parabolic and, as such, can be numerically integrated in a step-by-step 
procedure along the'boc 3 y surface. In order to cast the equations into 
a form, in which the step-by-step procedure can be efficiently utilized, 
the derivatives with respect to | and U are replaced by finite 
difference quotients. 

The method of linearization and solution used in the present 
analysis closely parallels that of Flugge-Lotz and Blottner (ref. 15) 
with modifications su gg ested by Davis and Flugge-Lotz (ref. I6) to 
improve the accuracy. These modifications involve the use of three- 
point implicit differences in the I -direction ■vjhich produce truncation 

* V , 

errors of order Alg) rather than (Ax) as in reference I5. The 

* V 

primaiy difference between the present '"development and that of refer- 
ence 16 is that the solution is obtained in the transformed plane for 
arbitrary grid-point spacing in' the Ir direction -and for a spacing in; 

i , ’ 

the T) -direction such that the spacing between ky two successive grid . 

t ' 

points is a constant. To the author's knowledge,; this numerical solu- 
tion technique has not been previously applied to either transitional 
or turbulent, boundary -layer flows. 
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The three -point implicit difference q.uotients are first developed 
for the particular grid-point spacing utilized herein. These results 
are then substituted into the governing equations for the conservation 
of momentum and energy in order to obtain a system of difference equa- 
tions. The most important, characteristic of the system of difference 
equations, other than. the fact that th^ are linear, is that the trans.- 
formed' normal component of velocity, V does- liot appear explicitly as 
an unknown at the solution station. Consequently,, the J 5 T -1 linear 
difference equations can he simultaneously solved to yield the N-1 
unknown values of 'F and 0 . Having obtained the H -1 values of F 
and 0, the TT-1 values of V can be obtained by a numerical integra- 
tion of the continuity equation. Some details of the implicit method 
have been purposely omitted; for these the reader is referred to 
Fliigge-Lotz and-Blottner (ref. 15). The reader' interested in a thorough 

t 

discussion of the various schemes that could be utilized to obtain either 
pure difference equations, as in the present paper, ‘ or’ difference- 

V 

* r ' ■* 

differential equations such as used. by Smith ‘and Clutter (refs. 8l and 

82) is referred to- Davis and Fliigge-Lotz (ref. I6). The- advantages and 

>' ' * 

disadvantages of inplicit differences in relation to explicit differ- 

i * ^ * • 

ences are discussed by Fliigge-Lotz and Blottner (ref. 15 )- and will not 


be discussed herein. 
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8.1 The Implicit Solution Technique 

8.1.1 Finite Difference Mesh Model 

It has heen shown for laminar houndaiy layers that equally 

i 

spaced grid points can he utilized in the normal coordinate direction 

(for example, see refs- l^^ and 15)*- However, for transitional and 

turbulent houndaiy layers, the use of equally spaced grid points is 

not practical for the present solution method. As previously mentioned 

in Section the viscous sublayer thickness is on the order of 

0.0016 to 0.015. 3h order to obtain a convergent (valid) solution to 

the governing equations, it is necessanry to have a number of grid 

points within the viscous sublayer. The viscous sublayer problem is 

t 

discussed in detail in Section where the errors resulting from 
improper grid-point spacing in the wall region. are demonstrated with 
a numerical example. Since a number of grid points -mast be positioned 
within this relatively thin layer in relation to the total thickness, 

6, one cannot afford- to utilize ‘equally spaced grid points from the ' 

I I * < 

viewpoint of computer storage requirements and processing time per test 
case. For example, if the total thickness of the T)-strip -v/as 100 for 

if ' ^ 

a particular problem -of interest, then the requirement of equally 
spaced grid points would mean on the order of 1,000 points in relation 
to 200 for variable spacing. In order to avoid the necessity of using 
an inefficient number of grid points, a variable grid-point scheme 
must be' utilized. 

The grid-point spacing in the Ti-direction used in the present 
analysis assumes that the ratio, of any two successive steps is a 
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constant j that is, the successive form a geometric progression. 

* * Y • » 

There are any number of schemes' that could be employed to distribute' ^ 

the points across the iq-strip) however, the pre.s'ent method has been 

^ ' 

found to be flexible and of utility. ' ‘ ' 

The desirability of having variable grid-point spacing in the 

■ ( I 

I -coordinate has previously been mentioned in Section 7.5* The' 
importance becomes clearly apparent for problems in which either the 

"i * » 

rate of change of the edge conditions is large or discontinuous bound- 
ary conditions occur, or in the -transition region where the mean pro- 
files are changing rapidly. A good example of the former case would 
be a slender, blunted cone in supersonic flow. Variable grid-point 
spacing in the I -direction could be utilized by having very small steps 
in the stagnation region where the pressure gradient is. severe 
(•favorable) and in some downstream region where transitional flow exists. 
A good example of the case of discontinuous boundary conditions would 
be a sharp-tipped cone with a porous insert at same downstream station 
through which a gas is being injected into the boundary layer. Rela- 
tively large step sizes could be utilized upstream of the ramp injec- 
tion; however, small steps must be used in the region of the ramp 
injection. Downstream of the porous region, as the flow relaxes, 
larger step sizes could be used. It is veiy important that small grid- 
point spacing be utilized in the transition region where the mean pro- 
files are a strong function of the intermittency distribution, lypical 
examples of the above-mentioned cases are shown schematically in 
figure 8 . Therefore, because of the requirements imposed "by a 
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completely general problem where the possibility exists for abrupt or 
rapidly changing edge conditions and boundary "values^ as well as the 
particular problem associated with the transitional and turbulent 
boundary layers.^ -variable grid-point spacing is utilized in the present 
technique in both th.e I and "Q directions. The grid-point spacing 
in the I -direction is completely arbitrary.. The grid-point spacing -in 
the T|-direction is such that the (i = 1,2, • ,PT) form a geo- 
metric progression. • . ' 

t 

In constructing the difference quotients, the sketch of the 
grid-point distribution presented in figure 9 useful for reference. 
The dependent variables S’ and © are .assumed known at each of the 

* t ' r 

N grid points along the m-1 and m stations , but unknown -at 
station m-i-1. The and Alg values., not specified to be equal 

are obtained from the specified x-values and equa- 

tion (7* 553')- The relationship between the Aflj^ for the chosen grid- 
point .spacing is gi-tf-en by the following- equation: 

■ ATI = (K)^"^ AT] (8.1) 

r X -1- 

In equation (8.1), K is the ratio of any two successive steps, At]^ 
is the spacing between the second grid point and the -wall (note that 
the first grid point is at the wall.), and DT denotes the total number 
of grid points across the chosen rj-strip. The geometric progression 
spacing of the grid points defined by equation (3*1) is but one of any 
numbe.r of possible schemes. This particular approach has been found to . 
be satisfactoiy by a number of investigators; for example, Bushnell and 



'6 



Terence grid, model. 
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Beckvdth (ref. 5 I) and Cebeci (ref. 78 ). The total thickness -of the 
■Q-strip can then be expressed as follows: 



At]_ 


1 - (K) 
1 - K 


N-1 


(K ^ 1) 


( 8 . 2 ) 


Another particularly important and useful relation is< 

V 



,(8.3) 


The spacing of the grid points >isv-coinplete2y 'determined for a specified 

{ ' ' * ' ^ 

set of and IT values. -The selection of ithe optimum .• 'K' and 

• - > ' ■ ' . • , 

N values for a specified depends upon the particular prpblem 

under consideration. The main objective in the sdlection is to obtain 
the minimum number of grid points with which a convergent solution may 
be obtained. That is^ in order to minimize the computer processing 
time per test case^ it is necessary to minimize the total number of 
grid points across the boundary layer and to maximize the spacing be 
between the solution stations along the surface. The laminar layer 
presents no problem since -a K value of unity is acceptabldj however^, 
for transitional and turbulent layers., the value of K will be a 
ntuaber slightly greater than unity, say l.C4. 

8.1.2 Finite Difference Relations 

Three-point implicit difference relations are used to reduce the 
transformed momentum and energy equations (eqs. (7.39) 9J^d (7.4o)) to 
finite difference form. As previously mentioned in Section 8.1.1, it 
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is assumed that all data are Imown at hhe stations m-1 and m (see 
fig. 9)- We then wish to obtain the unknown quantities at the grid 
points for the m+1 station. The notations G and H are 
utilized in the following development to represent any typical variable. 

Taylor series expansions are first written about- the unknown grid 
point (m+l,n) in the | -direction as follows: 


2 • 3 - 

•A^2/ - 

*^m,n“ ^itL+l^n " Wl,n 




and 


• f ' -(AI3. + ' 

*^m-lj,n “ '^m-M^n- “ (^^1 Wl,n **■ T Wl,n 


(Alp + AI2)' 
6 . 


^*^|||^m+ljn 


« • * 


(8.4b) 


where subscript notation has been utilized to denote differentiationj 

that is, G. = ( ^ 

S \.S5/ 

Equations (8.4a) and (8.4b) can be solved to yield 


/$g\ _ ^l^m+l,n ~ ^2*^m,n ^^^m-l,n 

\d|/m-i-l,n 2 A|g 

A|2(Ali +Al2)' . ^ 

^ ■* ' 


(8.5) 


and 


A|,AU/ A|o\ 

Sn+l,n = Vm,n ‘ 


+ 


(8.6) 
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The coefficients appearing in equations (8.5) and (8.6) 

are defined as follows: 


and 


Xi = 2 


+ 2 A^2 

f AI2 


<8.7) 


= 2 


All 


( 8 . 8 ) 


X^ = 2 


A|^A§2 


A5j_(A|i + A^g) 


(8.9) 


X), = 


_ ^^1 -*- ^2 


A|- 


( 8 . 10 ) 


5 A^ 


X. = 


( 8 . 11 ) 


Taylor series expansions are next written about the unknown grid 
point (m+l;n) in the 't]-direction as follows; 


A\ 


^+l,n+l “ Ai1n,(CJTi)m+l,n — 5— (G'T)T])m+l,n 


+ ^(g ) + 

6 ''^TiTiT] Wl^n 


(8.12a) 


and 



8o 


^ 2 
All T 

n-1/ 




flkriccj ) ■ ' + 

(■ ''^TjTiTi V+l,n * 


Equations (8.12a)’ and (8.12b) can' be solved to yield 


and 


^STi^/m+l^n. 


^l*^m+l,n+l " ^2^ai+l,n 


(*Vl - ^1n) 


TiTin 


'8g\ 

■^1 ^4^m+l,ntl ” " ^6*^m+l,n-l 


A% 


•^TiTiq + • • ' 


The . . . ;Yg coefficients appearing in. equations (8.13) 

are defined as follows: 


Y. = 


1 At)^(Ati^ h- At]^_3_) 


^2 = 


^^n-1 


Y^ = 


5 A%„i(A% + A%_q) 


^ (8.12b) 

I 

t 

(8.15) 

■ . (8'. Ill-) 
and (8.14) 

(8.15) 

(8.16) 


(8.17) 



8i 


AT) 


^1). = 


'n-1 




(8.18) 


and 


S' 


'^Vi - 
*1n ^Vl 


(8.19). 




’All, 


'n 


* ^Vl> 


( 8 . 20 ) 


For tlie case of equally spaced grid points in the’’ ^ '• and ‘ t) - 
coordinates^ equations (8.7) to (8'.^)-*and (8.15) to (8.20) reduce to 
the follomng relations': ' ' . 


X'2 ii; 

X5 .-I: 

\ =-'^ 
X5 = 1 


;(8.'21a) 


and 


Y =-^ 

1 AT)2 


Y2=2Yi 


Y3 = Y^ 


A 


X4 


(2 All) 


(8..21h) 
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•where A| and An represent the spacing het-ween the grid points in 
the I and n coordinates, respectively. 

Equations ( 8 . ' 5 ), ( 8 . 6 ), ( 8 .I 5 ), and ( 8 . ill-) can then he ifritten 
for constant grid-point spacing -as follows: 


,S|/m+l,n 


^^nna,n ~ ^^m,n ^m-l,n 
2 A| 







( 8 . 22 ) 


^m-i-l,n “ ^^m,n “ ^m-l,n ^ 


(8.25) 


^2 \ 
8 G ‘ \ 


^itt+l,n+l “ ^m+l,n Si+l,n-l At)^ 


^^8n^/m+l,n 


G, 


An" 


12 


+ • * • 


(8.24) 


and 


i'f) 

on/m+l,n 


*^m+l,n+l " ^ifl+l,n-l An^ ^ 
2 An " 6 


(8.25) 


Equations (8.22) to (8.25) are recognized -as the standard relations for 

equally spaced grid points. (See, for example, ref. I 6 .J 

• /^\ ■ ■ - ' 

Qjian-bities of the form 'that appear in the' governing 

■equations must he linearized in order to obtain a system of linear 
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difference equations. Quantities of this type are linearized "by 
utilizing equations ( 8 . 5 -)' and ( 8 . 6 )^ that is.^ 

\ °^/m+l,n 


+ 0(A^3_ 


(8.26) 


The procedure used to linearize nonlinear products such as 

/8g\/Sh\ 

(— l(— - j is the same as that used "by Flugge-Lotz and Blottner (ref. I 5 ) 

v^vv^v •. 

and is as follows: 


'SgV^h'^ 


\8T|/ IQt)] V 3i1/ ISt)]-'' : 

' 'm,n\ /m+l^n \ Ym,n\ /m,n 


V ^ V m, n\^ V m+1^ n 


. J ( 8 . 27 ) 


where the' terms- and are evaluated from egua- 

\OT]Jm,n \oT|ym,n, 

tion (8.14).;. but at the known station, m. Equating G to H in 
equation (8.27), the linearized form for quantities of the tsrpe 

^q,\2 

t:— is obtained? that is, - ■ 

on/ 


' mV 

M 




_ >/ 5g \ 

\SVm+l,n 


where 

\ovm+l,n 


m+l,n 

is obtained from equation (8.l4). 


(8.28) 


The preceding relations for the difference quotients produce 
linear difference equations when substituted into the governing 
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differentiai equations for the conservation of momentum (7»59) and 
energy (7.40). !Ehese linear difference equations may he written as 
. follows: 




+ Eln®m+l,n + ^^n^m-fl^n+l = 


^n^m+l^n-l ®^n^m+l,n ^^n®m+l,n- 




The coefficients Alj^^^Blj^, . . . . . . ,G-2^ are functions of quan- 

tities evaluated at stations m and m-1 and are therefore known. 

These coefficients are presented in the appendix. 

8.1.5 Solution of Difference Equations 

The proper boundary conditions to be used with the difference 
equations for the specific problem under consideration axe specified 
in equations (7*^5)- The system contains exactly 2(lir-l) mutualily 
dependent equations for 2 (n- 1) unlmowns since the boundary conditions 
are specified at the wall (i = 1 , see fig. 9)- This system of equations 
can then be simultaneously solved for the 2(lT-l) uhknowrs at the m+1 
station. The system is of a very special type;, since a large number 
of the coefficients in the system are zero.- The simultaneous- solution 
technique has been discussed in detail by ‘Fliigge-Lotz and Blottner 
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(ref. 15)- A discussion of the technique as applied to the incom- 
pressihle laminar boundary -layer equations is. also presented by 
Schlichting (ref. 2, pp. l8l-l8il-). '(See, also, Richtinyer, ref. 83.) 


8.1.4 Solution of Continuity Equation 

It is in^ortant to notice that the transformed normal component 
of velocity, V does- not appear as an unknown at station m +1 in the 
difference equations (eqs. (8.29) (8.3O')). This arises because of 

the way quantities of the tjrpe linearized (s.ee eq. (8.26)). 

Consequently, equation (7.38) can be numerically solved for the N -1 
unknown values of V at station m +1 once the values of F and © 
are known at station m+ 1 . Equation ( 7 * 38 ) can be integrated to yield 
the following relation for V at the grid point (-m+l,n): 


V 




m+l,n m+1,1 



(8.31) 


where represents the boundary condition at the wall and is 

defined in equation (7 '^ 7 ) as a function of the mass transfer at the 
wall, (p*v*)^. The integral appearing in equation (8.3I) can then be 
numerically integrated across the t]- strip to obtain the lT-1 values, 
of V. In the present analysis the trapezoidal rule of integration was 
utilized. However, any sufficiently ,accTArate numerical procedure could 
be used (see, for .example, refs. 84 or 85). 
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8.1.5 Initial Profiles 

Initial profiles for sta,rting the finite difference scheme’ are 
required at two stations since three-point , differences are utilized 
(see fig- 2}. The initial profiles at the stagnation point or line for , 
hlunt "bodies, or at x = 0 for sharp-ripped bodiesi are obtained "by an 
exact numerical solution of the ^similar boundary-layer equations. The 
equations- are solved by a fourth-order Adams-Bashforth-Moulton fixed- 
step size integration method with a fourth-order Eunge-Kutta technique 
(ref. 85 ) used to start the integration. The N-1 values of F, ©, 
and V which are now known at the W-l equally spaced grid points 
are numerically redistributed to U"-l grid points whose spacing is 
determined from equations ( 8 . 1 ) to (’ 8 . 5 ) if a variable spacing is 
required. The second initial -profile- located at station m is assumed 
to be identical to the one located at station m-1. .Any errors that 
might be incurred because of this assurption are .minimized, by using an 
extremely small AS, that is, an initial step size in the physical 
plane on the order of Ax = 1 X 10“^. The solution -at the unknown 
station, m+1, is then obtained by the finite difference method. One 
advantage of variable step .size in the I -coordinate is clearly demon- 
strated for hlunt body flows. Davis and Flugge-Lotz (ref. I 6 ) found 
that from the standpoint of accuracy and computer processing time 
step sizes on the order of O.OO 5 in the physical plane were required. 
However, in the present solution, extremely small, equally ^ spaced step 
sizes can be- utilized in the stagnation region and then increased to 
more realistic values once the errors due to the starting procedure 
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have approached zero. Davis and Flugge-Lotz (ref. l6) were primarily 
interested in obtaining solutions for only a few nose radii downstreamj 
consequently, equally spaced but sufficiently small grid-point spacing 
in the x-coordinate could be utilized. However, for engineering calcu- 
lations where solutions may be required as far downstream as 1,000 or 
more nose radii, one cannot afford to use equally spaced grid points. 

It is also advantageous to have the capability of starting the 
solution from experimentally measured profiles, especially in the case ■ 
of turbulent flow. Eiis capability has also been incorporated into 
the digital computer program used in the present analysis. This capa- 
bility is extremely useful for cases where one cannot easily locate the 
origin of the boundary layer, for example, nozzle walls. 

8.1.6 Evaluation of Wall Derivatives 

> ”” , * 

The shear stress and heat transfer at the wall (eqsp(7-97.) 

> « ■ ■ 
and (7.106)). are directly proportional to the' gradient of F and ^ 0 

. V 

evaluated at the wall, respectively. Using G to. represent a general 
■quantity, where Gjii+1,1 specified to be zero, the four-point ' 

difference scheme used to evaluate derivatives at the wall is as 
follows : 

(§)k+l,l " '^7^m+l;l + Wm+1,2 

where the coefficients , Y-j_q are defined by the following 

equations ; 
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^T=- 


(1 + K + k^)^[k(i k) - + (i + k) 

(1 + K)(l + K + k2)k5 At]]_ 


(8.53a) 


^ 8 ^ 


(1 + TC + K^) 

k 2 All-, 


(8.53b) 


and 


^9 


(1 + K + K^); 
(1 + KjK^All^ 


(8.53c) 


- 


(1 + K + K^,)k 5 AT]-, 


(8.53a.) 


For the case of equally spaced grid points in the i]-direction 
(K = 1),^ equations ( 8 . 53 )' become 


Y = --S— 
T (6 All) 

Yo - iQ 

8 (6 AT]) 


^9 r ('6 ATI) 

_ 2 


‘•10 


( •6 AT] ) 


(8.51^) 




and equation (8.52) reduces to the familiar four-point relation) that is, 


(ILa * - 


(8:. 35) 
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8.1.7 Eddy Viscosity Distribution 

The major advantage of the present solution technique in rela- 
tion to Cebeci, Smithy and Mosinbis (ref. 7^) ^md Beckwith and Bashnell 
(ref. 86), for example, is that the momentum and energy equations 
(eqs. (8.29) and -(8.50)) are simultaneously solved without iteration, 
whereas in the above two references the momentum and energy equations 
are each individually solved and iterated for convergence. The eddy 
viscosiiy function distribution 6 and e (see eqs. (7*15) (T*i6)) 

and their derivatives with respect to t] represent somewhat of a prob- 
lem unless extreme care is used in the initial extrapolation of the 
known values of %_i^n %,n unknown station, m+l,n. 

During the development of the digital computer program, the 
n\jmerical method would frequently become unstable in either the tran- 
sitional or turbulent flow region (see Section 8.1.8 for discussion of 
convergence and stability). This problem would always occur in one of 
two ways. In some instances an apparently converged solution would be. 
obtained, but the distribution of boundary-layer thickness would not 
be smooth. In other instances, ^diere the transition was abrupt or 
where boundary conditions were abruptly changed, the' solution would not 

converge. The problem -was finally traced back to its origin, which was 

* 

a nonsmooth or "rippled" eddy viscosity distribution across the layer. 

I , , ' 

These "ripples" first occurred in the region where, the inner and outer- 

eddy viscosity models were matched. If the initial "ripples" were 

* ’ I ^ * 

below a certain level, the solution would . apparently converge, but 

* I , . • 

slightly nonsmooth boundary -layer thickness. distributions would occur. 
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If the initial "ripples" were above a certain level, the disturhance 
would grow very rapidly as a function of | and propagate throughout 
the layer as the solution proceeded downstream. VJhen this occurred, 
no valid solution -could he obtained downstream "of the initial 
disturhance. 

ihe initial extrapolation of the known values of 
^m-l,n unknown grid point (m+l,n) is obtained as follows (see 

eq. (-8.6)): 

, t • 

■ Wl,n “ ^5^m-l,.n , , '(-8-56) 

However, there is no assurance that the distribution of the extrapolated 
values at. station m+1 will be smooth across the- layer for a.l 1 possible 
flow conditions. If "ripples’" occur in the extrapolated £ distribu- 

< I 

tion and if these "ripples" are of sufficient magnitude to cause the 

I 

— * J ■ 

sign of the derivative of e with respect to > to alternate, then 

t 

the method becomes highly unstable. , 

The requirement of small grid-point spacing in the law of the 
wall region contributes to the instability problem in that the size of 
an "acceptable ripple" is a function of the g-rid-point spacing being 
utilized in the t] - direction. For turbulent layers where the viscous 
sublayer is. relatively thick, the grid-point spacing in the outer 
portion of the law of the wall region will be large in comparison to • 
cases where the viscous sublayer is relatively thin. Consequently, 
the former case can tolerate a "larger ripple" in the region of the 
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match point than can the latter case -without experiencing an alterna- 
tion in the sign of the derivatiTe of "e with respect to f]. 

There are t-wo possible -ways to eliminate the problem caused by 
the "ripples" in the eddy viscosity distribution. The first approach 
is to develop an iteration scheme in -which case the present solution 
technique has no advantage in relation to the technique used in 
reference 86j that is, the advantage of the unique simultaneous solu- 
tion would be lost. The second approach is to numerically smooth' the 
extrapolated eddy viscosity distribution prior to the matrix solution. 
Both approaches were tried by the author during the development phase 
of the digital computer program. !Hie second approach was incorporated 
into the solution and will be discussed in the remaining portion of 
this section. 

The problem posed by the "ripples" in the eddy -viscosity dis- 
tribution, -if they exist, can be avoided by utilizing a three-point 
mean value for e at station m+l,n, that is. 


C^av)m+l,n “ 


%+l,n-l + %+l,n- + ^+l,n+l 


(8.57) 


where denotes- the three-point mean of the eddy viscosi-fy func- 

tion. In the present analysis the eddy viscosity functions appearing 

V 

on the' right-hand side of equation (8*37) first obtained at each’’ • 
grid point across the m +1 station from equation '(8.36). . Having 
obtained these values, the three-point mean is -e-valuated at each of 
the H-l grid points from equation ( 8 . 37 '). The matrix solution is 
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then obtained for equations (8.29) (8-30)' Having obtained the 

W-1 values for Fj ®, and V at station m+1, the eddy viscosity 
distribution is recalculated at the m+1 station prior to moving to 
the next ^ grid-point station. This procedure has been found to be 
stable under all circumstances and to yield convergent solutions for! 
transitional and fully turbulent boundary layers.' 

8.1-. 8 Convergence and Stability 

The itiglicit difference scheme utilized- in the present, analysis^ 
is .consistent^ that is, it may be considered as a fomal approximation 
to the governing partial differential equations. A finite difference 
scheme is consistent if the difference between the partial -differential 
equation and the difference equation goes to zero as the grid-point 
spacing approaches .zero. That is, a scheme is consistent if the 
truncation error goes to zero as the grid-point spacing approaches 
zero. The difference quotients used in the present analysis satisfy 
this requirement. However, the fact that a system is consistent does 
not necessarily imply that the solution of the difference equations as 
the grid-point spacing approaches zero converges to the actual solution 
of the partial differential equations. In' order for convergence to be 
assured, the difference equation system must be convergent. Unfortun- 
ately, for the boundary-layer equations as used herein there is no 
■completely satisfactory mathematical analysis -with which the conver- 
gence or divergence of the present scheme can be demonstrated (ref. 15) 
The convergence of the system can best be s'tudied numerically by vaiy- 
ing the grid-point spacing. This approach is used in Chapter IS of 
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the present paper. Convergence prohlems cati become very critical in ' 

the turbiilent region of the boundai:y layer when the viscous sublayer is 

( ■ 

not properly treated- Of course, one should keep in mind that there 
are no exact solutions of the mean turbulent boundary -layer equations 
against which comparisons for convergence may be made. 

The stabili'hy of the finite difference scheme, as opposed to the 
convergence, can be carefully studied mathematically. The three -point 
implicit scheme used in the present method is stable regardless of 
the grid-point spacing, (ref- 15^ PP- The stability of the 

method is the major advantage of the implicit, difference scheme in 
relation to the explicit approach -where it is necessary to control the 
grid-point spacing through a stability parameter. (See, for example, 
refs. 87 to 9I-) 

8.1.9 Vis CPUS' Sublayer - 

"The relative thinness of the viscous sublayer {.Section 7»^*l) 
and its importance in determining the correct shear stress and heat 
fl"ux at the -wall requires that very small grid-point spacing in the 
T] -direction be ’utilized in the wall region for the present n-umerical 
procedure. Although a careful study of either the- minimum number 'Of 
grid points that must be located in the viscous sublayer or the optimum 
spacing of the, grid points across the entire boundary lay.er has not 
been .completed, it has been found by the author that at least 5 to* 10 
grid points in the- viscous sublayer -at each x-solution station will 
yield a convergent solution. (This study is currently being completed 
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and will le reported in a subsequent UASA publication containing the 
digital computer program. ) 

Preliminary data obtained from experimental studies and a 

thorough literature survey by the Gas Dynamics Section at the Langley 

.1. 

Research Center indicates that the viscous sublayer thickness, y , 

s. I • 

ranges from approximately 8 to' l 4 over the Mach number range from 3 

to 47 for adiabatic to cold wall conditions and for air or helium flow 

mediums. (These data are currently being, processed- for publication 

« * 

W. P. Harvey and F. L. Clark at the Langley' Research Center.') Conse- 

H. 

quently, a mean value 2 ~ equation' (7.116a). allows a 


reasonable festimate of y* , 

. s. Z. 


as follows.; 
t 


V* = -- ^ 


S • Z « 


T* 

_Si 

p* ’ 

w 


. (8.38) 


i , 


where the local values of and in equation (7jH6a) have been 

replaced with the kno'wu wall values in equation (.8.38).' An accuxaie 
estimate of ' the value' of t* for a particular case c^ be obtained 
from a ■mmiber of current -correlations, for example, Raiding and Chi 

(ref. 52). Therefore, an estimate of y* , can be obtained, from 

s. i . 

equation {'8.38)] however, the numerical solution is obtained in the 
transformed plane. The .sublayer thickness in the transformed plane can 

then be obtained directly from equations (7.53t>) and the estimated y* 

s * ^ « 

value. The final step is to estimate the boundary-layer thickness and 
the optimum total number of grid points across the boundary layer. 
Having, reasonably accurate estimates of the viscous sublayer thickness, 
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the total hoxmdary -layer thickness^ and the optimum nimiber of grid 
points^ equations (8.1) to (8'. 5) can then be used to obtain the spacing 
parameter^ K, In order to place from 5 to 10 grid points in the vis- 
cous sublayer., An example of the effect 'of grid-point spacing on the 
convergence of the numerical method is presented in Section 9*1* 

8. 2 Digital Computer Program. 

The repetitious, nature, of the implicit finite-difference scheme 
makes, the method well suited for digital computers. The present solu- 
tion technique i/as programed by the author for the CDC (Control Data 
Corporation) 6000 series computer operating on Scope '5. The coded 
program will not he presented herein since it is quite lengthy as well 
as currently being- documented for publication by the national Aeronautics 
and Space Administration. 

The main objectives in the development of the computer program 
used in the present study were accuracy^ minimum computation process 
time, flexibility, and simplicity. The accuracy requirement is, of 
course., of paramount importance in any numerical procedure. - The 
CDC 6000 series machines cost on the order of $1 thousand per hour to 

* i 

operate. Consequently, the conputer process time per test case must 

be held to a minimum "by optimum programing procedures. Both of these 

requirements should be carefully considered in any digital computer 

program. However, flexibility and simplicity are-' often overlooked when 

the primary objective of the study is to obtain solutions for one 

» • 

particular class of flows in order to satisfy one specific requirement. 




The ir^iortance of simplicity in input requirements' is also often- neg- 
lected when the needs of other personnel associated with the research 
center^ university or aerospace industry are hot taken into eonsidera- 

t I * 

tion- These objectives were necessary in the present case' since the 

( f 

program was developed by the author for -general use at the national' 

Aeronautics and Space Administration, Langley Research Center. Some of 

the more ingjortant points concerning the flexibility '6f‘ the program 

will be discussed in the remaining portion of this section. 

The input requirements that must be specified in order to obtain 

the extern^ flow 'are discussed in Section 7.8. The input values for 

I* ^ CO-* ^ 00 ^ ^ dependent upon only 

the free-stream test conditions and flow medium. However, V* which 

e 

must be specified as a function of the boundary-layer coordinate, x*, 
is a function not only of the test conditions and flow medium, but also 
body geometry. The pressTire di-stribution, P|, as a function of the 
boundary -layer coordinate, x*, is obtained from experimental data when 
possible or from the method presented in reference 80, when applicable. 
A very useful option for engineering studies of the boundary -layer 
characteristics has been provided by programing Modified Newtonian 
pressure relations (ref. 95) directly into the program for blunt-body 
flows; that is, ■ . 



(8.59) 
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where the Biaxitiium. pressure coefficient^ 



may he escpressed as 




7 + 5 


,ax 


7 + 1 


1 - 


!^(.7 + 3). 


(8.4o) 


■niis option is particularly useful for cases where engineering design 
data are required for. a particular class of configurations over a 
hroad range of flight conditions. A good example would he the class 
of power-law bodies for \diich equation (8.59) will yield an ^accurate 
representation of the true pressure distribution at hypersonic .speeds. 

1 ^ s ? 

(ref. ^). This class of bodies may be described mathematically as 
follows : ' ' ■ 



‘(8.ia) 


i , , 

where r^ designates the base radius for a body of length L.% For ‘ 
this particular example, the power-law exponent, a,' could be varied 
to study the viscous drag over a range of flight conditions with a’ 
mlnimTora. of input requirements. Another very useful application would 
be to study the heat transfer to the vehicle over a given flight trajec 
tory. Di this case the variable input required for the solution would 
be specific values of ’^oo points along the' tra- 

jectory. The most probable location of transition 'and the extent of 
the transition region worild be obtained from empirical correlation 
relations (see Section 7.5). 

It should be noted that the above-mentioned class of bodies 
(eq. ( 8 .i^l)) presents somewhat of a problem in that the boundary -layer 
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coordinate, x, cannot be obtained in closed form as a function of the 
body coordinates, z and r^^. This presents no problem in the 
present program, since the program contains an iteration subroutine 
which will obtain x as a function of the body coordinates for any 
geometry where r^ can be expressed as a function of z. 



IX. EXMPtE SOLUTIONS FOR THE SYSTEM OF EQUATIONS 

The finite difference -sol-ution technique developed, in Chapter VIII 
together with the transition criteria^ transitional flow structure, and 
transition extent discussed in Chapter VII is applied In the present 
chapter to a number of typical supersonic- and hypersonic configurations 
of current interest. In cases the gas is assumed to be perfect air 
with a constant ratio of specific heats, 7, equal to 1.1)-, a constant 
Prandtl number, o, equal to O.72, and a constant static turbulent 
Prandtl number, equal to 0,9. The molecular viscosity, q , is 

evaluated from Sutherland's viscosity law (eq. (7.32)). The external 
pressure distributions used are either experimental or were obtained by 
the author from an- exact .solution of the full inviscid- Euler equations 
obtained at the Langley Research Center (ref. 80). 

The selection of a typical set of test cases always presents a 
problem since there are. so many possibilities from which to. choose. 
However, the cases considered in the present chapter should serve as an 
indication of the merits of the solution technique as well as the 
validity of the system of equations. 

9.1 Hi^ Reynolds Number Turbulent Flow 

The accurate prediction of boundary-layer characteristics for 
high Reynolds number turbulent flow is., io^ortant in the design of high- 
speed vehicles.- In particular, it is important- to be. able to predict 
■with accuracy the skin friction drag. An excellent comparison case of 
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high Reynolds number turbulent flow is the data of Moore and Harkness 
(ref. 95)* experimental skin friction data were measured with a 

floating element type balance. The test model was a 10-foot-long, 
ij— foot-wide sharp (0.002 to 0.004 inch thick) leading-edge flat plate. 
The test conditions were as follows': 

.K> = 2-8 

P* = 2.088 X 10^ Ib/ft^ 

t,,oo 

^t,~ = 5*6 X 10^ °R 

~ = 9-^7 X 10 "^ 

The. experimental transition location was not reported in 
reference 95- Consequently, for the numerical calculations, the tran- 
sition location was determined by the stability index (Section 7'5*2.l) 
and was assumed to occur at the x- station ’where achieved a value 

of 2V500. The extent of the transition region, was automatically calcu- 
lated from equation (7.85). The intermlttency distribution was calcu- 
lated from equation, (7.81). The solution was started at the leading 
edge of the sharp flat plate (x =' 0) by obtaining ^ .'exact 'numerical 
solution of the similar boundary -layer equations (eqs. (A-47) to (A-^9)) 
by the method discussed in Section 8.’l-5* The grid^point spacing was 
varied in both the I and t] directions in order to check for conver- 
gence. It should be noted that the selection of (^iax^cr “ 2,500 

' - ' ■ ‘V * 

was based partially on experience and parti^ly ph the results presented 
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in Section S-Q -where the variation of (^ax^cr 1!^* , is dis- 

cussed. !Qie location of transition ranged' from ,0.8ij- to 1.4 inches 

’ . ' ’ . * - 

from the ,leading edge for (^ax^cr tjet-ween 2,500 and' 3 ^ 500 , respec- 
tively. This variation in transition location had negligible effect on 

i * ' • 

the skin friction coefficients' in -the> fully -turbulent region, of flow. 


The numerical results for the skin frictidn. coefficient distri- 
bution are con^ared to the experimental ^ da-fa a in .figure 10 (a). The 
agreement is excellent over bhe entire Reynold's n-umber range of the 
experimental data for K = 1.04; however, for K = 1.01, convergence 
was not attained. It should be noted at this point that the terms 


convergence and stability, as used in relation to the numerical me-fahod, 
are defined herein as in reference 15 and are discussed in Section 8.1.8. 
The divergence (fail-ure to obtain a convergent solution) for K = 1.01 
is attributed to an insufficient number of grid points in the wall 


region; in particular, in the thin viscous- sublayer region. This 
particular problem is discussed in .Section 8- 1.9* The effect of the 
grid-point spacing parameter K on -the numerical solution -was studied 
for additional values of 1.02, 1.05, 1.05, and 1.06. The solution was 


found to diverge for K < 1.02 and converge for K > 1.02. The 
results, for K > 1.02 :( other than 'K =■ 1.04,)- are not presented since 
the variation of -with K would not be discernible "if plotted to 

the scale of figure 10 (,a-). The convergence criteria used in this 
parti c-ular. example was to decrease the grid-point spacing until any 
change which occurred in at a given x-solution station was beyond 

the fourth significant digit. The laminar curve shown in figure 10(a) 
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■was obtained by suppressing transition (r = 0). Grid-point spacing 
In the I -direction was varied from Ax = 0.001 to 0.,04> and convergence 
to the accuracy of the experimental data was obtained for all values. 
Because of the abruptness of transition, step, .sizes greater than 
Z&: = 0.04- were not studied. 

A comparison to an experimental veloci'ty profile is presented 
in figure 10 (b ) . ‘ The profile -was obtained at an Ee^j. value of 
8.4-5 io"^. No experimental data wera avallablp for values less 

P' ' ' * i 

than 10 „ The agreement between the numerical results and the experi- 
mental profile is seen to be v6iy good. " ^ 

, ; . 4 » ■ 

Typical profiles of E, ©, prior 

to transition (Re„ = 1.4-6 X 10^), mid'way through’ the transition region 
(Re^ = 1.96 ,X 10^), and at the completion of transition (Rex = 5-44 X 10^) 

i 

are presented in- figures 10(c) to 10(f), respeetivelyl ITo^ experimental 
data were available in this region^ however, the results are, in 
general, "bypical of any transition region (see Section 9*S). The 
intermittency distribution and bomdaiy -layer tbicRness distribution 
in the transition region are presented in figures 10(g) and 10(h), 
respectively. 
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(e) Transitiou region total pressijre profiles 
Fig-ure 10.- Continued. 



(f) Transition region Mach number profiles 


Figure 10.- Continued 
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9*2 Tripped Tuicbulent Boimdary Layers 

In many of the current supersonic mnd-tunnel facilities it is 
often necessary to trip the laminar houndary layer artificially in 
order to obtain turbulent data. A good example of turbulent data 
obtained in tripped boundary layers -is that of Coles (ref. $ 6 )'. These 

» . I 

data .were obtained in the Jet !^opulsion Laboratory 's _20-inch super- 

sonic wind tunnel. The test model, was a sharp leading-edge flat plate. 

/ > 
r 

The free-stream Mach number was yaried from;1.966 to. Test 

numbers JO, 20, and 22 (see. page JJ of ref. 96) were- selected ‘as typical 
comparison cases. For these three eases the laminar boundary./ layer was 
tripped by a fence located at .the leading edge -of the flat, .plate (see 
fig. Jl-Q of ref. 96). The skin friction was measured at three surface 
locations with a floating element balance. Boundary -layer, profiles 

' ■ t 

were measured at x* = 21.1)-8 Inches. 

The test conditions for the three comparison cases are listed 
in tabular form helow. 


Coles test No. 

nni 




50 

1.982 


ESI 

0.8295 

20 

3.701 

2868'. 5 

. 561 

.7152 

22 

■4.55^ 

8132.2 

554 

.6764 


Transition was assumed to occur near the leading edge of the 
plate for the numerical calculations,, x.^. ^ 


= ,0.005 ft, and to be 
















io8 


completed (r = 1.0) at „ =? 0.010 ft. Twenty equally spaced 
x-solut±on stations were used in the region 0 < x* < 0.010 ft^ the 
x-solution stations were then equally spaced 0.01 ft apart over the 
remainder of the plate. The total number of grid points in the 
T)-direction and the grid-point spacing parameter K were assigned 
values of 50I l-Cit-, respectively. The -computer processing time 
per test case was approximately 3 minutes. The numerical results are 
compared to the experimental velocity and Mach number profile data in 
figures 11(a) to 11(f) for the three test cases. The agreement between 
' the numerical results and the experimental data is very good for all 
- three test cases. In particular, it should be noted that the experi- 
mental skin'friction coefficients (see figs. 11(a), ll(c), and 11(e)) 

■ were predicted to within 1 percent, -which is well within the accuracy 

f 

range of 2' percent as quoted for the data in reference 96. 

Humerical results obtained hy Pletdier^ (ref. 97J see, also, 
ref. 98) are also presented on figures ll(a) to U(f). A three-layer 
mixing length concept was used in reference 97 "to model the turbulent 
transport terms. The resulting equations were solved by an explicit 
finite difference, technique using the DuPort-Frankel (ref. 99) differ- 
ence scheme. The agreement between the numerical results of refer- 
ence 97 -and those of the present solution is very good for the three 
test cases. To the author’s knowledge, the method used in reference 97 
is the only explicit finite difference technique currently being used 
to solve the mean flow turbulent boundary -layer equations. 

^The author would like to thank Professor R. H. Fletcher, 
Department of Mechanical Engineering, and Engineering Research 
Institute, Iowa State University, Ames, Iowa, for permission to use 
these data. 
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(a) Velocity profile and skin friction coefficient for = I. 982 . 

Figure 11.- Comparisons -with expe 3 :djaental data for tripped turbulent 

boundary layers. 



(b) Mach number profile for kto = I. 982 . 
Figure 11.- Continued. 


(c)' Velocity profile and skin friction coefficient for 

1^ = 5.701. 

♦ 

figure 11.- Continued. 
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9*5 Laminar Flow With Mass Injection 

In some instances it may tecome necessary to protect high- 
performance vehicles by transpiration cooling {mass infection at the 
wall). An example of laminar data that may he used for comparison is 
presented hy Marvin and Akin (ref. 100). The data were obtained over 
a range of injection rates for a sharp-tipped 5° cone. The cone i^as 
solid for X* < 3*75 inches j the remaining portion of ■ the cone was 
porous. The test conditions were as follows; 


= 7-h 


Pf = 8.0J- X 10^ Ib/ft^ 

t,“ ' 


Tf = 1.5 X lo5 
t^w 





= 3.8 X 10"1 


The, air ‘injection rate ranged from I.O 56 x 10“^ slugs/ (ft^-sec) to a 
maxifflum value of 3 * 8 ^ x 10 “^ slugs/ (ft^-sec). 

A comparison of the heating rate at the wall normalized by the 
heating rate at the wall just prior to the ranp injection (x = l) is 
presented in figure 12(a). The mass injection distribution was not 
uniform, as can be seen from the plot of the actual injection distribu- 
tion normalized iy the ideal injection rate. The non-'uniform distri- 
bution of mass injection was utilized in the numerical solutions. The 
agreement between the numerical results and the experimental heat trans- 
fer data is excellent over the entire injection range. 
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(Clw) 




X 

(a) Comparisons to heat-transfer data. 

Figure 12.- Hypersonic laminar boundary layer flow 
with mass injection. 


h 


For the maximum injection rate ((p*v*).jj. = 5 >8^ slugs/ft^/sec);, 
the boundaiy layer separated from the cone surface at x = 1 .. 92 * ®ie 
calculated velocity and Mach number profiles at x = 1-92 are presented 
for the four test cases in figures 12(b) and 12(c), respectively. 


Present solutions , slugs/^ 

0 

LOSexio-"^ 

2.112 I 

3.840 ^ 


x= 1.92 



(b) Velocity profiles. 
Figure 12.- Continued. 
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(c) Mach number profiles. 
Figure 12. - Continued. 


The vorticity Reynolds number distribution for (p*v*)^ =0 is pre- 
sented in figure 12(d). The maximum value of the vorticity Reynolds 
number^ \ietx^ seen to Increase with increasing x. In particular^ 
the y/s value at which occurred is in excellent agreement with 

the location of the critical layer position (see ref. 6j). The effect 
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of mass injection on the vorticity Reynolds number distribution is 

presented in figure 12(e). Increasing mass injection is seen to 

increase the value of x-station as well as’ move the 

location at which the maxinrum occurs j (y/6)v toward the outer 

^ax 

edge of the boundaiy layer. To the author's knowledge no data are 
currently available for the critical layer position with mass injec- 
tion at the wall. 

This particular example is a case where variable step size must 
be utilized in the region of the ramp injection (x = 1.0). For this 
test case, the Ax values (grid-point spacing in x-direction) were 
constant and equal to 0.01 up to x = 0.99 s-i' which point the step 
size was decreased to a value of 0. 001 through the ramp injection 
region. The step size was then progressively increased up to a value 
of O.Ol at X = 1 . 5 . The flow was laminar; consequently, K was set 
to unity and 101 equally spaced grid points were utilized in the 
■n-direction. The digital computer processing time per test case was 
approximately 2 minutes. 

It should be noted that the test conditions as listed in 
references 100 and 101 are not correct. The conditions given herein 
are correct and were obtained by the author from Mr. Joseph Marvin at 


the Ames Research Center. 




(e) Vorticity Reynolds number distribution for variable - 

mass injection. 


Figure 12.- Concluded. 


119 


9*^ Blunt Body Flo-w 

The numerical procedure as presented in the present paper has 
"been applied to a number of blunt body supersonic and hypersonic flows. 
A -typical case is presented in figure I3 for hypersonic flow. The 
pressure distribution was obtained by the author at the Langley Research 
Center (ref. 80). The test conditions were as follows: 

= lO.ij- 

= 2.2l|-9 X io5 ib/ft^ 

.Solutions for this particular example have been obtained by Marvin and 
Sheaffer (ref. 101) and Clutter and Smith (ref. 102). The distributions 
of the heating rate at the wall referenced to the value at the stagna- 
tion point (x = 0) and the- shear stress at the wall are compared to 
the resTxlts presented in references 101 and 102 and are presented in 
figure 13. 

Equaldy spaced .grid points were used in the I -direction. Solu- 
tions for Ax values of 0.001 to 0.02 were obtained to check for con- 
vergence. Oscillations in the due tothe starting procedure 

approached zero after approximately 10 to 15 steps along the surface. 

The data shown in figure I3 were obtained using Ax'= O.OO5, although 
much larger steps could be used for engineering calculations., say. 

Ax = 0.01. The grid-point spac^g in the -q-direction -was constant • 

\ ’ 

(K = 1) with a value of 0.01. The computation time per body station.' 

' i j 

‘ 1 

was. about 3 seconds as compared to- 42 seconds .in. reference 102. 


T*^„= 2.2 X 10^ °R 
— = 2.5 X 10"^ 
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Figure I 3 .- I^ersonic Hunt Taody flow. 


9*5 Highly Honsimilar Flow With Transverse Curvature 

Turbulent boundary layer data for flows where variable ^pressure 

* ■ - ' 

gradients exist are few in number. 'A good example of a case where 
both favorable and adverse pressure gradients occur ‘as well as where 
transverse curvature effects are important \is‘ the data of Winter., 

r S ' - I ‘ 

Rotta, and Smith (ref. IO 5 ). The model used in the study was an' axi- 

t ’ 

symmetric piece-wise continuous, configuration .and is presented in 
figure 1 ^ 1 - (a). . ' ; ■ . ‘ 

Experimental data are presented in referenc.e IO 3 only for the 
region where the boundary layer is turbulent; hojwever, the -solutions 
presented herein -were obtained by starting the calculations at the tip 
of the sharp cone fbrebody '(x = 0-). This particular configuration has 
received considerable attention over the past '2-year period. 




z/c' = 0; c*= 5 ft 

Tq/c = 0.36397 (z/c) 

z/c = 0,14158 

Tq/c =6.3473 [z/cf- 6.492EI (z/cf+1.22668 (z/cf+ 0.33498 (z /c) - 0.00461 
z/c = 0.45725 ■ 'r 

Tq/c = - 0.236382 lz/cf+ 0.60762 (z/cf- 0.59552.1 (z/c)+0.22694 

z/c = 0.64473 . 

r^/c = - 53.58258 (z/c) + I50.87806'(z/c) -158;03866(z/c) +72.96859(z/c)- 
• - - -12.4939 

z/c?0.763i7 

Tq/c = 0.236382(z/cf+a396789/z*/c)^+0,290242(z/c)-0.06400l 
z/c = 1.0 


(a) Geometry of configuration. 

Figure l4.-- Comparisons with data for highly nojasimilar supersonic 
flow with transverse curvature effects. 
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CalctCLations have been 3nade hy Herrihg and Mellor (ref. 10^)-), Ceheci, 
Smith, and Mosinskis (ref. 7^), and Bushnell and Bec3cwith ’ (ref . 51 )• 
However, the above-mentioned solutions were a ll started by utilizing 

, f ’ . 

the experimentally measured profiles at station z'= 2. To the author’s 

^ ^ 1 

knowledge the solutions presented in the- present paper are the first to 
be obtained without any dependence whatsoever- on experimental profile 
or skin friction data. ^ ‘ 

The test conditions for, tjie two cases considered are ast follows: 
Case 1 • Case 2 ' 

= 1-598' H».= l.70 


Pf =■ 9. 21 X 10^ Ib/ft^ 
Tf = 5-56 X IG^ °R 
= 9-76 X 10”^ 

i| ►X* 


P* = 9-92-x 10^ Ib/ft^ 
t,«> 

T*" = 5.56 X 10^ °R 

t,<» ■ 


rri*- 

t,<» 


= 9.71 X 10 


-1 


The experimental Mach number distributions are presented in' figure ll(b). 
The edge Mach number distribution was used as an input to the digital 
computer program instead of P*. Equations (7.I51) to (7.15I) and 
(7-31) were used to calculate the req}iired edge conditions from the 
specified edge Mach number distribution. The initial conditions behind 
the conical shock wave were obtained from reference 105. 

The initial profiles required to start the finite difference 
solution were obtained 'by an exact solution of the similar boundaiy- 
layer equations (eqs. (A-47) "to (A-l9)) at x = 0. Transition was 
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initiated at the solution station, ■where Xj^ax achieved a value 

of 2,500* The transition extent -was then automatically computed from 
equation (j.6j). The grid-point spacing in the | -direction varied 
from a maximum Ax value of 0.01 to a minimum value of 0.001 in the* 
regions of large pressure gradients. Variable grid-point spacing in 
the Ti-direction -was required with a K value of 1.04. Calculations 
were made for K values of I.05 and I.05 to insure convergence. The 
computer processing time per test case -was approximately 4 minutes. 

The numerical results are con^ared with the experimental data 
for momentum thickness and skin friction coefficient distributions in 
figures l4(cj to l4(f). 

The agreement, between the numerical -and experimental momentum 
thickness and skin friction coefficient distributions is very good for 
both test cases. In particular, note the agreement -with the minimum 
Cfg data point in transition (fig. l4(f)) and the erratic heha-vior 
of the data in the region x < 1 (figs. 14(e) and 14(f)). It is also 
of interest, to note that while the solutions with transverse curvatuxe 
were in closest agreement -with the 0 values for x < 3*5^ the solu- 
tions -without transverse curvature agreed best with the values 

in the same region.. A similar trend is presented -in reference 78. 











0 


(e) Skin friction coefficient for ly^ 
Figure l4.- Continued. 
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9*6 Adverse Pressure Gradient Turtulent Flow- 

The design engineer is often required to make estimates of the 
skin friction, heating rate, and other boundary-layer characteristics 
for turbulent boundary -layer flows in. adverse pressure gradient regions. 

Good examples of such flows would be the boundary layer over a deflected 
control surface (conpression surface) and over the centerbody of a 
supersonic engine inlet. An example of turbulent data obtained in an. 
adverse pressure gradient is presented by McLafferty and Barber 
(ref. 106). One of the test configurations was a flat plate with a 
single-radius curvature compression surface at the rear of the plate. 

The test conditions were as follows: 

Mco = 5-0’ T* = 6.10 X 10^ °R 

t,«> 

' rr^ 

P* = 2.116 X 10^ Ib/ft^ = 9.28 X 10“^ 

t,“ . rrtJt V 

For this particular exauple the pressure gradient across the 

i 

boundary layer was significant. However, for the' present- numerical 

•i 

calculations, this factor was neglected since the^system. of- equations' 
as used herein does not contain the y-momentum equation. The experi- 
mental Mach number at the edge of the boundary lajrer was used ‘to ' obtain 
the required edge conditions. This -distribution is presented in 
figure 15(a). The station x = 0 marks the beginning of the com- 
pression surface. No experimental: measurements- of the skin friction 
were made. 

For the present solution the calculated momentum thickness 
Reynolds number at the beginning of the compression, x = 0, was 




X 


(a) Experimental Mach number distribution. 

Figure 15.- Comparisons to experimental data for adverse pressure 
gradient "supersonic turbulent flow. 
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, . i 

matched to the experimental value of Ecq = 2,5^0. Equally spaced 

’ 4 

x-solution stations -were used; Ax = 0.001. The total humher of grid 
points in the ii-direction and' the grid-point spacing parameter K were-_ 
assigned values of 201 and 1.04-^ respectively. The computer processing 
time -was approximately 2 minutes. For this particular test case it was 
necessary to increase the value of Kg from 0.0168 to O-O672- in the 
outer region eddy viscosity model. This increase in Kg was necessary 
in order to obtain attached flow throughout the compression. It should 
he noted that this particular example is the only case in the present, 
paper where- the value of Kg is different from tliat given in 
Section J.h.l.; that is. Kg = 0. OI68. The numerical- solution resvilts 
for momentum thidsness and form factor are presented in figure 15(h). 

The calculated momentum thickness distrihution agrees well with the 
experimental data for -2 < x < 0. 5; however, the numerical resTlLts 
fall somewhat helow the data throughout the remaining compression. The 
agreement between the solution and the experimental form factor distri- 
hution is good throughout the coupression. 

The eddy viscosity model as used herein (see Section 7-^) 
developed from data obtained in zero pressure gradient flows. The 
effect of streanline curvature has been studied by Bradshaw (ref. IO7). 
Bradshaw found that streamline curvature (centrifugal effects) had a 
strong effect on the turbulence structure for flow geometries -tdiere the 
ratio of the boundary -layer thickness to the radius of curvature of the 
boc3y surface was on the order of I/500. For concave- streamlines, such 




X 


(Id) Comparisons to eDcperimental momentum thickness and fonn 

factor distributions. 


Figure 15 •- Concluded, 
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as in the present case;, it was shown (ref. IO 7 ) that the mixing length 
should he increased in order to account for the effect of circvature on 
the turbulence structure. 

This particular’ test case has also been studied by Beckwith 
(ref. l4).. Beckwith avoided separation in his solutions by computing 
the velocity gradient, dUg/dx, from the static pressure distribution 
at the edge of the boundary layer, P^, rather than’ P^, as used in 
the present solution; however, u^ and Pg were calculated from P^ 
as in the present solution. It should be recalled that separation, was 
avoided in the present solution by increasing the eddy viscosity in 
the outer region. It was s-uggested in reference 14 that the premature 
separation problem should be corrected by the inclusion of the 
y-momentum equation and that more data should be obtained in adverse 
pressure gradient flows in order to develop appropriate mixing length 
■or eddy •viscosity models. 

9.-7 Flow Over Sharp -Tipped Cones 

Fisher (ref. 56 ,) studied the effect of the unit Reynolds nmber, 
E*, on transition for a 10° half -angle cone at a free-stream Mach 
number of 7 . These data pro-vide an excellent case for comparison as 
well as to exhibit the flexibility of the digital computer program. 

The test conditions were as follows: 
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^=0 = 7 ^ 

2-89 X 10^ to 8.7T X 10 ^ it/ft^ 

T* = 1.05 X 105 °R 

t,“ 

= 5.2 X 10 

mX- 

t,“ 

The ‘botindaiy-layer edge values were obtained from, reference 105* The 
•experimental location of transition and the extent of transition -were 
used in the numerical calculations. The conputer processing time per 
test case was approximately 2 minutes. 

Comparisons of the numerical results to the experimental , 

Stanton nimiber distributions are presented in figures l6(a) through 
16(f). The value of (^^max^cr noted on each figure. The agreement 
between the numerical results and the es^erimental data is very good 
for a.T 1 six test cases. In particular^ note that the' numerical .solu- 
tion predicted the correct peak heating and overshoot characteristics 
of the t-wo' hipest unit Reynolds number cases ,(;see figs. l6(-e) and 

16(f)): , ;i 

Two cases similar to the above ones are presented in figures l6(g) 
■and 16(h). These data were obtained by P* Cal-vln, Stainback (unpublished) 
at the Langley Research Center on a 10° sharp cone model. -The test 
conditions for these data were as follows: 




* 6 

(a)' " 5.00 X 10' . 

Comparisons with experimental Stanton number distributions 
for hypersonic flow over sharp-tipped cones. 



© Exp. data 
. Present solution 

(>^max)cr= ^450 



(h) 1^, = 7; rJ = 5 . 9 ^ X 10^. 

Figure 16 .- Continued. 






6 8 10 12 
X* in 


(c) = 7, rJ = 5-30 X 10^. ' . ■ . 

Figure l6.- Continued. 



2 4 6 8 10 12 


X , in. 

(d) Ho = 7. Rg = 6.69 X 10^. 
Figure l6.- Continued. 
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Case 1 

Moo = 8 

P* = 2.91 X iq5 

T* = 1 . 458 . X 1q 5 °R 
t,“ 

^ = 4.0 X io"l 

These two cases are typical of approximately 15 cases calculated liy 

I ► 

the author from the same source. Results from these remaining ‘cases 

will be presented in a, future publication. They are presented herein 

* 

because of the length of turbulent flow in relation to that obtained 
by Fisher (ref. 56 ). The (^ax^cr Is indicated on the .figures. 

The agreement in the peaic heating region is excellent; in particTzlar, 
note the overshoot characteristic and its agreement ‘with the data. 

9 . 8 Flow Over Planar Surfaces 

O’Donnel (ref. IO 8 ) studied laminar, transitional, and turbulent 
boundary-layer flows over a hollow cylinder. Velocity profiles were 
measured at various stations 'along the cylinder. The test conditions’ 


Case 2 


= 8 


P^ - 5-65 X 10^ Ib/ft^ 

T^ ■ = 1.367 X lo5 9 r 

t,“ ' . 


3 ^ 

= 4.2 X lO”! ‘ 

in* \ - 


< ‘ 


Hx, = 2.41 

'Pf = 4.95 X 10^ to 8.49 X 1q 5 Ib/ft^ 


were as follows: 




For this particular set of 'calculations, the experimental tran- 
sition location was utilized’; however, the extent of transition was 
calculated from equation ,(7.85). Consequently, the only variable 
inputs to the computer program were the specific values of the total 
pressure, P* and the transition location, 

The velocity profile comparisons are presented in figures 17(a) 

6 

through 17(e). For a unit Reynolds number of 0. 672X10 , the boundary layer 
was laminar throughout the measured area. The velocily profiles are 

similar and the agreement between the numerical results and experimental 

data is very good (see fig. 17(a)). For” a unit ^Reynolds number of 

« 

/C ■' r ‘ ’ 1 , • 

2.88 X 10°, laminar, transitional, and turbulent flow occurred (.see'* ' 

* I ^ 

fig. 17(1)). Similar results were obtained for a unit Reynolds number 
of 5.76 X 10^ (see fig. 17(c)). For unit Reynolds numbers of 8.64 x10^ 
and 11.5 ,X 10°, the flow was turbulent as- presented in figure l?(d.) and ‘ 

i 

^ 4 • 

17(e), respectively. The profiles are seen to be similar for the tur- *- , 
bulent region. Contparisons with the experimental momentum thicknesses 
are presented in figure 17 (f). The agreement is seen to be excellent 
over the entire unit Reynolds number range i 

A further example of laminar and transitional boundary-layer 
flow for hypersonic test' conditions is presented in figure 18. 



m 



(a) = 0.672 X.IO^. 


> , 

Figure 17*- Comparisons with, velocity profile data and momentum 
thickness for laminar, transitional, and turbulent flow over 
a hollow cylinder at =’2.4l.' ... 



(b) rJ = 2.88 X 10^. 


Figure 17-- Continued- 



Ik2‘ 



Pigiire 17.- Continued. 



(d) H*,= 8.64 X io^.. 


Figure- Continued.^ 

} 





These data were obtained by Johnson (ref. 109) -on a shai^ leading -edge 
flat plate model. The experimental test conditions, were ^as follows: 

= T.8 

p* =■ 8.827 X 10^ ib/ft^ 

Tf = 1.45 X 10^ °R 

^ = 5.88 X 10".l 
t,“ 

The agreement between the mmerical results and the experimental 

Stanton number distriboition is "very good. 

The e3q)erimental transition location, x* . , was- used in the 

t^x 

calculation; however, the transition extent was' calculated from- equa- 
tion (7.87.). (note that the transition from laminar to turbulent flow 

was' not coir5)leted; that is, x* „ ■was not reached. ). The total number 

, t,r 

of grid points in the T)-direction and the .grid-point spacing parameter, 
K, were assigned values of 201 and i.04, respectively. The con5)uter 
processing time was approximately 2 minutes. 

This particular test case is presented to emphasize that the 
stability index, (\iax)cr ^ invariant as suggested by Rouse 

(ref. 64) but is, in fact, a strong function of the unit Reynolds 
number, R*. For a.Tl of the previous test cases, (^ax)'cr 
varied over a relatively small range 2,500 < (')^ax)cr 
which might lead one to assume that, while not invariant, (^ax^cr 
varies only slightly. However, the present test case, value of 4,000 
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considera63y 'extends this range from 2^500 to. 4,000. The ntmerical 
results obtained "fc^r assuming that transition would occur at the 
x~station where equal to 2,500 is shown on figure l8. This 

I, 

calculation, clearly indicates that no assumptions concerning the loca- 
tion of -experimental transition should be made on the basis of the 
magnitude of the stability index as being an invariant or a near 
constant. 

The stability index is presented as a function of unit Reynolds 
number for a number of experimental studies in figure 19. These values 
were generated by obtaining the boundary -layer solutions corresponding 
to the experimental' test conditions and plotting the value of 
at the experimental transition location as a function of the unit 
Reynolds number of the experiment. Some of these cases have been pre- 
sented as test cases in the present paper (refs. 109, 5^,?, 108, and 
two Stainback cases). A number of points of interest -concerning the 
stability index are apparent in figure 19 • For example, these data can 
be represented by the' general equation 

(^^ax’^cr ” ^4 ^°Sl0 ^e ^5 (9*l) 


where and represent functional relations of the transition 

4 5 

parameter phase space. It should he noted at this point that all of 
the data presented on figure 19 represent similar boundary -layer flows 
and that nonsimilar effects such as pressure gradients mi^t consider- 


ably alter these trends. 




Figure 















The limited experimental ■ data presented on figure ^19 in^cate 
that is invariant, at least, for similar flows since the slopes of 

the curves for the data from references 56, IO8, and the data from 
Stainhack are identical. However, the coefficient 0 .^ is definitely 
not an Invariant for these data. The data from reference 56 and 
Stainhack (unpublished) were obtained at nearly equal edge and wall 
boundary conditions; however, the value of is considerably 

different. This suggests that is a strong function of the test 

environment since the data were obta,ined in different wind-tunnel test 
facilities; that is, perhaps the acoustic environment would appear in 
The data presented for the adiabatic cylinder (ref. IO8) exhibits 
the same dependence on as that of reference 56 and Stainhack 

(unpublished); however, since these data were, obtained for adiabatic 
wall conditions and at a much lower Mach, number than that for the 10 ° 
sharp cone, would appear to be a function of at least Mg, T^/Tg, 

■and the test environment. The test environment would include the 
effect of acoustic radiation, free-stream turbulence, and temperature 
spottiness of the flow. 

The author is currently studying transition data obtained over 
a broad range of test conditions for both similar and highly nonsimilar 
flows. The basic approach being followed is simple in concept; that is, 
obtain highly accurate boundary -layer solutions using the present solu- 
tion technique for the test conditions under which the experimental 

) 

data were obtained. However, while simple in concept, the process is 
quite lengthy since care must be exercised in the selection pf the 



test cases. It is important that only well-documented data he studied 
in which the local test conditions as well as the test environment are 
carefully reported. This represents a problem in itself since in, much 
of the currently existing transition data no attempt was made by the 
experimenter to measure the basic characteristics of the test facility 
in which the data were obtained (test environment). The eventual goal 
of this study is to obtain not only more transition correlations but 
also more insight into the actual transition process. 



X- Disajssiosr.AM) conclusions 


A system of equations which describe the laminar, transitional, 
and turbulent compressible boundary layers for either planar or axi- 
symmetric flows as well as a numerical method by which the system can 
be accurately . solved has been presented. 

The transitional boundary layer structure was developed from the 
point of view of- the statistical production and gro-vfth of turbulent 
spots. The passage of these spots over points on the surface results 
in an alteration of laminar and fully developed turbulent flow. -These 
alternations are described by an intermittency factor which represents 
the fraction of time any -point spends in turbulent flow^ or the prob- 
abili-by at .any given instant of time that a specific point will be 
engulfed in a turbulent spot. The intermittency factor was used to 
modify the turbulent transport models developed for fully turbulent 
flow. 

Stability theoiy and its relation to experimental transition has 
been briefly reviewed. Experimental transition data were discussed 
and methods were suggested that should allow reasonable estimates to 
be made for the most probable location of transition and the extent of 
the transitional flow region. The primary objective of the present 
paper -was to present a system of equations ,and a solution technique 

V , 

with which the boundary -layer characteristics could be obtained regard- 
less of -whether the flow was laminar, transitional, ‘or turbulent. ^ 
Consequently, it -was necessary to discuss, in' some detail, the 
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transition, problem since the location- of transition for a particular 
flow geometry must ’be either specified from experimental data or calcu- 
lated from some correlation equation. A conplete treatment of the 
transition problem would be a study in itself that would require many 
years to complete. The treatment of the transitional flow structure, 
as used in the analysis, is not exact in the sense of obtaining a sys- 
tem of time-dependent equations from- the full Navier-Stokes. system, 
but is instead based upon existing experimental data. A thorough treat- 
ment of the transition location and the transitional flow structure 
still remains as one of the major unsolved problem -areas in fluid 
mechanics. However, the method as presented herein does provide 
accurate predictions of the mean flow profiles and wall values of heat 
transfer and .skin friction in the transition region. To the best -of 
the author's knowledge, this represents the first time this procedure 
has been incorporated into a solution technique for the complete non- 
similar boundary-layer equations. ' 

The turbulent boundary, layer was treated by a two-layer concept ' 
with appropriate eddy viscosity models being developed for each l^er 
to replace the Reynolds stress terms in the mean turbulent boundary— ■ 
layer equations. A specifiable static turbulent Prandtl number ’was 
used to replace the turbulent heat flux term. A constant turhuleht 
static Prandtl number was utilized. However), any model ‘could he 
directly employed once developed from experimental data. 

The numerical method used to solve the generalized .system of 
equations is a three-point implicit difference scheme for variable 



152 


grid-point spacing in "both spatial coordinates. The method is self- 
startingj that is, it requires no experimental data input. The method 
has been found to be highly efficient with regard to flexibility, 
digital computer processing time, and accuracy. The method is 
inherently stablej that is, no constraint is placed on the grid-point 
spacing by -a step size stability parameter such as in the case of 
explicit finite difference schemes. However-, the grid-point ■ spacing 
is restricted by the maximum truncation error that c^ be accepted for 
a given test case. To the best of the author’s knowledge, this repre- 
sents the first time' this particular solution technique’ has been 
applied to transitional and t'urbulent boundary -layer flows. . ' 

A number, of test cases have been presented and compared with 
■experimental data for supersonic and hypersohic flows over planar and ■ 
axisymmetric geometries. These test cases have included laminar, 
transitional, and turbulent boundaiy-layer flows ‘with both favorable 
and mild adverse pressure gradient histories as' well as a case of 
highly adverse pressure gradient flow. Mass injection at the wall and 
transverse curvature effects were also considered. The agreement, 
between the numerical resiilts and the experimental data was very good 
for all of the test cases eonsidered with the exeeption of the ease for 
high adverse pressure gradient turbulent flow. The agreement for this 
particular case was fair; however, much more work needs, to be done in 
this particular area. In particular, for adverse pressure gradient 
flows approaching separation, the normal momentum equation must be 
included in the governing system of equations’. There is also a current 
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need for well-documented experimental data for transitional and turbu- 
lent boundary-layer flows with, pressure gradient histories. Hiese data 
are needed to develop mixing length and eddy viscosity models for 
pressure gradient flows. One of the more important parameters that 
should be measured is the intermittency as a function of both spatial 
coordinates. One of the main advantages of the present numerical 
method is the conceptual sin^licity that can be utilized in the formu- 
lation of the turbulent flux terns. (This advantage would be true for 
all Implicit finite difference techniques. ) This advantage allows the 
use of "numerical experimentation" through which mixing length and 
eddy viscosity' models may be efficiently studied in relation to experi- 

r 

mental data. Through this "numerical experimentation" procedure the 
knowledge, of basic, turbulent mechanisms could irltimately be improved. 
However,^ well-documented data must be available for this .process to be 
successful. , ' ' ■ 

■ . In concluding, it .has been shown that the system of equations 
presented herein and the. solution technique by which they are -solved 
can be used to obtain accurate solutions for laminar, transitional, 
and turbulent compressible bouhdary-layer flows. 
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XIII. APPENDIX 


■Equations (8-29) and ( 8 . 5 O) are the difference equations used to 
represent the partial differential equations for the conservation of 
momentuni ■ and energy; respective!!^. These equations are repeated for 
convenience as follows: 

■^n^m+l^n-l ®\L^m+l;n ^^n?nri-l;n+l ®^n^+l,n-l 

+ ,'Eln®m+l;n + ^’ln®m+l;n+l == (A-1) 

■^n^m-+l;n^l *^^n^Bi+l;n+l ®^n®m+l;n-l 

' *. ' 

^2n®m+l;n -*• == ^2^ (^-2) 

These equations are obtained from equations ( 7 . 39) and (7.4-0) and the 
difference quotients presented in Section 8.1.2. The coefficients -Aln.^ 
Bl^; etc.; in equations (A-l) and (A-2) are functions of quantities 
evaluated at stations m and m-1 (see fig. 9 ) and' are therefore 
known. TSiese coefficients are as follows: 



(A-3) , ■ 

L •* » 

•Kin = ^ 1^1 ~ Y 2 H 5 > f 'H 5 

(A-4) 

^ Cl„ - %% + XllSll ■■■ ■■ , 

' (A-3) 
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Din = "Wy 

y 

El„ = -5 D1 + 

n n D 

= - ^ “In 

= -2Wy 

Y. 

^6 

i 

\ 

^6 

1 

D2 = Y,K „ - YyH^ 
n 3 10 6 12 

/ 

E2j^ = - Y2 H^q - Y^H^2 

= %0 ^ ^4^2 


(A-6) 

(A-7) 

(A-8) 

(A-9) 

(A-10) 

(A-U) 

\ 

(A-12) 

(A’-13) 

\ 

(A-14) 


(A-15) 



l6j 

and 

G2^ = + Hg(FY)^ + Hg(T^)^ (A-i6) 

The coefficients . . . :,Yg, X^;,...;,X^, etc.j are functions of 

the grid-point spacing and are defined in equations (8.15) to (8.20) 
and (8.7) to (8.11), respectively. The coefficients etc., 

are defined as follows: 


n - . F (FT) 

1 m-fl ml 


= V ^ - L _ (E ^C’_ + E’ C ) 
2 ml ml 33 q 1 ml ml ml 


“3 = - Vml^ml 


Hi, =H3 


\n 


H = B F 
5 m+l ml 


H^ = -P 

6 m+l 


H = V - L ^(E ,C’ + E*.C _ ) 

7 ml ml ml ml ml mil 


(Arl7) 

(a-i8) 

(A-19) 

t 

(A-20) 

(A-21) 

(A-22) 


(A-23) 
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Hq ■ “<x , L E » C 
o m+1 ml ml ml 



-E L' -C ^ 
ml ml ml 


TJ _ TJ ^ 

^10 "^9 77 " 

■^ml 


Hu = 


and 


The imdefined quantities appearing in equations (A-I7) to 
are defined as follows: 


'ml 


= X|,F 
4 m^n 




m-l,n 


T 


'ml 


X, 0 - x^e ^ 

4 m,n 5 m-l^n 


V T -= X,,V - X^V , 
ml 4 m,n 5 m-l^n 


F = X F - X F 
m2 2 m,n ^5 m-l^n 


(A-2i|- ) 

(A-25) 

(A-26) 

(A-27) 

(A-28) 
(A-28 ) 

(A-29) 

.(A-30) 

(A-31) 


(A-32) 
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^ “ \/^ 



(Air only) 


L’_ 

ml 



(Air ©a^) 


^ml " ^^avWl,n Section 8.1. 7, eq. (-8.57))* 



V+l,n 


% 


^4^m,n+l 


- „ 

5 m;,n 


- T 

6 m,a-l 


(See -eq. '(■8.li)-,)) 




W^n+1 " ^5^m,n " ^6^m;,n-l 


(A-33) 

(a-34) 

(A-35') 

1 

(A-36) 

' \ 

(a-37)^ 

(A-38') 

(A-39) 

■(A-l^o) 


T_=:Yi0 -Y® - Y^0 

1 4 m,n+l 5 6 m,n-l 


(A-41) 
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and 


^m+1 


El 

u„ d| / . 

, e ^/m+1 


(See eqs. ( 7 .^ 4 -!)) 


-(A-42) 


V+l 



(A- 45 ) 


The transverse curvattire terms are contained in the quantities 
Cml and C^-, which appear explicitly in the H , H » H ^ Hg, and 
coefficients. The transverse curvature term in the transformed 
plane (see eq.. (7*4:5)) may be written as follows: 


t^^ = 1 + 


2oij (W) cos 0 


Pe^e 



(A-44) 


where t represents the ratio t/tq and is a known quantity for the 
N-1 grid points at station m-1 and m. Then, the extrapolated values 
at m+l,n are obtained as follows where the parameter C is used to 
represent t^^ : 


C = X, C - X^C T 
ml 4 m,n 5 m-l,n 


(A- 45 ) 


C’, = Y, C - y^C- - Y^C T _ 

ml 4 m,n+l 5 6 m,n-l , 

, i 


■ (A-46) 


Two quantities (symbols) as of now remain undefined., 'These are 
FT and W which appear in. equations (A-I 7 ) and (A-44), respectively. 
These are code symbols used in the digital program. 'The code symbol W 



171 


appearing in equation (a- 44) is used’ to either retain or neglect the 
transverse curvature terms for axisyrametric flows; that is, W = 1 or 0, 
respectively. For planar flows,, the transverse curvature term does not 
appear since j equals 0. 

The code symbol FT (flow "type) appearing in equation (A-17) is 
used to either retain or neglect the nonsirailar teims in 'the .governing 
differential equations; that is, FT' = 1 or 0, respecti-yely. I-f' FT 
is assigned a value of ’uni’by, the solution to the nonsimilar equations . 
(eqs. (7.58) to (‘7-40')) is obtained. If FT is /assigned a valiie of 
zero, the locally similar solution is obtained; that is, the following ' 
system of equations are solved. 

Continuity 

:~ + F = 0 ^ (A-47) 

dT) ' ■ ' , . 

. 1 V 

t 

Momentum ; ' - 

^ I - I) ° 

Fnergy 

V zi efef = 0 (Ajf9) 

• 5t;\ -a 8ti/ \^t\J 

The governing equations for the locally similar system are 
obtained from equations (7*5^) to (7*40) by neglecting derivatives of 
the dependent variables F, ©, and Y -with respect to the strearawise 
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coordinate, The capahility of obtaining locally similar solutions 

is desirable in that for a given test case the locally similar and 
complete nonsimilar solutions can be obtained for the identical pro- 
gram inputs and numerical procedures. Consequently, the effects of 
the nonsimilar terms on the boundary layer characteristics can be 
determined by direct comparison of the results. 



